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1. Introduction 

Among jet obscrvablcs, the three-jet production rate in electron-positron annihilation plays 
a very prominent role. The initial experimental observation of three-jet events at PE- 
TRA [1], in agreement with the theoretical prediction [2], provided first evidence for the 
gluon, and thus strong support for the theory of Quantum Chromodynamics (QCD) [3]. 
Subsequently the three-jet rate and related event shape observables were used for the pre- 
cise determination of the QCD coupling constant ag (see [4] for a review). Especially at 
LEP, three-jet observables were measured to a very high precision and the error on the ex- 
traction of Ug from these data is dominated by the uncertainty inherent in the theoretical 
description of the jet observables. This description is at present based on a next-to-leading 
order (NLO) calculation [5-10], combined with next-to-leading logarithmic (NLL) resum- 
mation [11,12] and inclusion of power corrections [13]. The calculation of next-to-next-to- 
leading order (NNLO), i.e. C(af), corrections to the three-jet rate in e+e" annihilation 
has therefore been high on the list of priorities for a long time [14]. 

Besides its phenomenological importance, the three-jet rate has also served as a theo- 
retical testing ground for the development of new techniques for higher order calculations 
in QCD: both the subtraction [5, 15, 16] and the phase-space slicing [7] methods for the 
extraction of infrared singularities from NLO real radiation processes were developed in 
the context of the first three-jet calculations. The systematic formulation of phase-space 
slicing [9] as well as the dipole subtraction [10] method were also first demonstrated for 
three-jet observables, before being applied to other processes. 
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Over the past years, many of the ingredients necessary for NNLO calculations of jet 
observables have become available: two-loop corrections to all phenomenologically relevant 
massless 2 — 2 [17] and 1 — 3 [18, 19] reactions were computed already several years ago, 
while one-loop 2 — 3 [20] and 1 — 4 [21] matrix elements are available for even longer. 
Despite all ingredients being available in principle, until recently it was not possible to 
perform NNLO calculations of any kind of exclusive observables, since techniques for the 
extraction of multiple real radiation singularities at NNLO were not sufficiently developed. 

Up to now, the only general method for handling this problem was the sector de- 
composition technique [22, 23] for the treatment of real radiation singularities. Using this, 
NNLO calculations of exclusive processes were performed for e+e~ — 2j [24], Higgs pro- 
duction [25] and vector boson production [26] at hadron colUders, as well as the QED 
corrections to muon decay [27]. 

Furthermore, exploiting the specific kinematic features of the observable under con- 
sideration, exclusive NNLO results were derived for the forward-backward asymmetry in 
e+e~ annihilation [28], for e'^e~ 2j [29,30] and most recently for Higgs production at 
hadron colliders [31]. 

In the present paper, we employ a recently developed general technique for the treat- 
ment of infrared singularities, antenna subtraction [32], to derive the NNLO corrections to 
three-jet production in electron-positron annihilation. The first phenomenological applica- 
tions of our results to the thrust distribution were documented earlier in [33] . 

The paper is structured as follows: in section 2, we outline the perturbative calculation 
of jet observables and summarise the antenna subtraction method used here. The imple- 
mentation of this method requires phase space mappings, which are described in Section 3. 
All relevant tree-level, one-loop and two-loop matrix elements arc listed in section 4. Sec- 
tion 5 briefly summarises how the NLO corrections arc implemented using antenna sub- 
traction. Sections 6-12 contain the subtraction terms constructed for all colour factors 
relevant in this calculation. The numerical implementation of all terms into a parton-level 
event generator is described in Section 13. As a first example of the implementation, we 
discuss the NNLO corrections to the thrust distribution in section 14. A summary and an 
outlook on applications is given in Section 15. 

2. Perturbative calculation of jet observables in e+e~-annihilation 

To obtain the perturbative corrections to a jet observable at a given order, all partonic mul- 
tiplicity channels contributing to that order have to be summed. In general, each partonic 
channel contains both ultraviolet and infrared (soft and collinear) singularities. The ul- 
traviolet poles are removed by renormalisation, however for suitably inclusive observables, 
the soft and collinear infrared poles cancel among each other when all partonic channels 
are summed over [34]. 

While infrared singularities from purely virtual corrections are obtained immediately 
after integration over the loop momenta, their extraction is more involved for real emission 
(or mixed real-virtual) contributions. Here, the infrared singularities only become explicit 
after integrating the real radiation matrix elements over the phase space appropriate to the 
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jet observable under consideration. In general, this integration involves the (often iterative) 
definition of the jet observable, such that an analytic integration is not feasible (and also 
not appropriate). Instead, one would like to have a flexible method that can be easily 
adapted to different jet observables or jet definitions. 

Three types of approaches for this task have been developed so far. Phase-space slicing 
techniques [7,9,35] decompose the final state phase space into resolved regions, which are 
integrated numerically and unresolved regions, which are integrated analytically. The sec- 
tor decomposition approach [22,23] divides the integration region into sectors containing a 
single type of singularity each. Subsequently, the phase space integration is expanded into 
distributions. In this approach, the coefficients of all infrared divergent terms, as well as 
the finite remainder, can be computed numerically. Finally, subtraction methods [5,15,16] 
extract infrared singularities of the real radiation contributions using infrared subtraction 
terms. These terms are constructed such that they approximate the full real radiation ma- 
trix elements in all singular limits while still being simple enough to integrate analytically. 

To specify the notation, we define the tree-level n-parton contribution to the m-jet 
cross section (for tree level cross sections n = m; we leave n 7^ m for later reference) in d 
dimensions by, 

da^ = A/- ^ d$n(pi, ...,Pn;q)^ \Mn{pi, ...,Pn)\^ J^\pi, • • • ,Pn). (2-1) 
n 

the normalisation factor J\f includes all QCD-independent factors as well as the dependence 
on the renormaliscd QCD coupling constant a^, denotes the sum over all configurations 
with n partons, d<I>„ is the phase space for an n-parton final state with total four-momentum 
in d = 4: — 2e space-time dimensions, 

(n) 

while Sn is a symmetry factor for identical partons in the final state. The jet function Jm 

(n) 

defines the procedure for building m jets out of n partons. The main property of Jm is 
that the jet observable defined above is collinear and infrared safe as explained in [10,38]. 
In general Jm^ contains and (^-functions. Jm'' can also represent the definition of the 
n-parton contribution to an event shape observable related to m-jet final states. 

|A^„P denotes a squared, colour ordered tree-level n-parton matrix element. Contri- 
butions to the squared matrix element which are subleading in the number of colours can 
equally be treated in the same context, noting that these subleading terms yield configura- 
tions where a certain number of essentially non-interacting particles are emitted between 
a pair of hard radiators. By carrying out the colour algebra, it becomes evident that non- 
ordered gluon emission inside a colour-ordered system is equivalent to photon emission off 
the outside legs of the system [36,37]. For simplicity, these subleading colour contributions 
are also denoted as squared matrix elements |A^mP, although they often correspond purely 
to interference terms between different amplitudes. The precise definition depends on the 
number and types of particles involved in the process. However, all colour orderings are 
summed over in ^^^^ with the appropriate colour weighting. 
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From (2.1), one obtains the leading order approximation to the m-jet cross section by 
integration over the appropriate phase space. 



d(TLO 



da' 



(2.3) 



Depending on the jet function used, this cross section can still be differential in certain 
kinematic quantities. 

2.1 NLO antenna subtraction 

At NLO, we consider the following m-jet cross section, 



da 



V 

NLO 



(2.4) 



The cross section da^j^Q has the same expression as the Born cross section da^ (2.1) above 
except that n m + I, while da^^Q is the one- loop virtual correction to the m-parton 
Born cross section do"^. The cross section daf^^Q is a (preferably local) counter-term for 
da^j^Q. It has the same unintegrated singular behaviour as da^j^Q in all appropriate limits. 
Their difference is free of divergences and can be integrated over the (m -|- l)-parton phase 
space numerically. The subtraction term dcr^^Q has to be integrated analytically over all 
singular regions of the {m + l)-parton phase space. The resulting cross section added to 
the virtual contribution yields an infrared finite result. 

Several methods for constructing NLO subtraction terms systematically were proposed 
in the literature [10, 15, 16,39-41]. For some of these methods, extension to NNLO was 
discussed [32,42,43] and partly worked out. Up to now, the only method worked out in full 
detail to NNLO is antenna subtraction [39,40]. In our calculation of NNLO corrections to 
thrcc-jct obscrvablcs, wc used this method, which we briefly outline in the following. The 
details of the method, and a full definition of the notation, can be found in [32]. 

The basic idea of the antenna subtraction approach is to construct the subtraction 
terms from antenna functions which encapsulate all singular limits due to the emission 
of unresolved partons between two colour-connected hard partons. This construction ex- 
ploits the universal factorisation of both phase space and squared matrix elements in all 
unresolved limits. The full antenna subtraction term is then constructed by summing 
products of antenna functions with reduced matrix elements over all possible unresolved 
configurations. 

At NLO, the antenna subtraction term thus reads: 

X ^ \Mm{Pl, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l)|^ JmHPl^ ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l) ■ (2-5) 



The key ingredient is the phase space mapping which relates the original momenta 
PiiPo-:Pk describing the two hard radiator partons i and k and the emitted parton j to a 
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Figure 1: Illustration of NLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + l)-particle phase space. The term in square brackets repre- 
sents both the antenna function Xfj^ and the antenna phase space d^Xi^t • 

redefined on-shell set pi,PK which are linear combinations of Pi,Pj,Pk- The phase space 
mapping yielding this redefinition is described in detail in section 3 below. With this 
mapping,the phase space factorises, 



d$m+l(pi, • • • ,Pm+i;q) = d$m(pi, • • • ,PI,PK, ■ ■ ■ ,Pm+l\q) ' d<^Xi^k{pi,Pj,Pk;Pl + Pk) ■ 



The other elements of the subtraction term also depend on either the original mo- 
menta Pi,Pj,Pk or the redefined on-shell momenta pi,PK but not both. This enables the 
subtraction term to completely factorise. 

To be more specific, both the m-parton amplitude and the jet function depend 
only on pi, . . . ,pi,PK, ■ ■ ■ ,Pm+i i-e. on the redefined on-shell momenta pi,PK- On the 
other hand, the tree-level three-parton antenna function X^jj^ depends only on Pi,Pj,Pk- 
Xfjj^ describes all of the configurations (for this colour ordered amplitude) where parton 
j is unresolved. It can be obtained from appropriately normalised trcc-lcvcl thrcc-parton 
squared matrix elements. The antenna factorisation of squared matrix clement and phase 
space can be illustrated pictorially, as displayed in Figure 1. Together particles i and k 
form a colour connected hard antenna that radiates particle j. In doing so, the momenta 
of the radiators change to form particles / and K. The type of particle may also change. 

One can therefore carry out the integration over the antenna phase space appropriate 
to Pi, pj and pk analytically, exploiting the factorisation of the phase space of eq. (2.6). 
The NLO antenna phase space d^Xi^^ proportional to the three-particle phase space, as 
can be seen by using m = 2 in the above formula. For the analytic integration, we can use 
(2.6) to rewrite each of the subtraction terms in the form. 



where |7W,rxP, Jm and d$m depend only on;ji, , . . . ,pi,PK, ■ ■ ■ ,Pm+i and d^Xijk 
depend only on Pi,Pj,Pk- This integration is performed analytically in d dimensions, yield- 
ing the integrated three-parton antenna function ^ijk, to make the infrared singularities 
explicit and added directly to the one-loop m-particle contributions. 



(2.6) 




d^X,,, Xf^, = \Mmf d^m Xijk 
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Figure 2: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles are 
colour connected. 



2.2 NNLO antenna subtraction 

At NNLO, the m-jet production is induced by final states containing up to (m + 2) partons, 
including the one-loop virtual corrections to (m + l)-parton final states. As at NLO, one 
has to introduce subtraction terms for the (m + 1)- and (m + 2)-parton contributions. 
Schematically the NNLO m-jet cross section reads, 



S 

NNLO 



+ 
+ 



ys,i 

NNLO 



da 



V,2 

NNLO ' 



(2.7) 



where daf^j^j^Q denotes the real radiation subtraction term coinciding with the (m + 2)- 
parton tree level cross section da^j^^^Q in all singular limits. Likewise, da^^]^Q is the 
one-loop virtual subtraction term coinciding with the one-loop (m-|- l)-parton cross section 
^^NNLO '^^ singular limits. Finally, the two-loop correction to the m-parton cross section 
is denoted by duj^l^j^Q. 

At NNLO, individual antenna functions are obtained from normalised four-parton tree- 
level and three-parton one-loop matrix elements. The full antenna subtraction term is then 
constructed by summing products of antenna functions with reduced matrix elements over 
all possible unresolved configurations. 

In dcr^^^Q, we have to distinguish four different types of unresolved configurations: 

(a) One unresolved parton but the experimental observable selects only m jets; 

(b) Two colour-connected unresolved partons (colour-connected); 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour- unconnected) ; 
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(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

Among those, configuration (a) is properly accounted for by a single tree-level three-parton 
antenna function like used already at NLO. Configuration (b) requires a tree- level four- 
parton antenna function (two unresolved partons emitted between a pair of hard partons) 
as shown in Figure 2, while (c) and (d) are accounted for by products of two tree-level 
three-parton antenna functions. The subtraction terms for these configurations read: 



^2 ^"^+2 



XI -^ijk \^m+\{Pl, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+2) 



'Jm^^Kpl^ ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+2) 



^^NNLO = X ^^rn+2{pi, ■ ■ ■ ,Pm+2;q)-^ 



m+2 



m+2 



y^ijkl ~ ^ijk^IKl - -^jkl^iJL) 

jk 

X \Mm{Pl, ■ ■ ■ ,PI,PL, ■ ■ ■ ,Pm+2)\'^ Jm\pi^ ■ ■ ■ ,PI,PL, ■ ■ ■ ,Pm+2) 
d'^NNLO = X '^^m+2{pi, ■ ■ ■ ,Pm+2; Q)-^ 



m+2 



'm+2 



X ^ijk ^mlK \Mm{Pl, ■ ■ ■ ,PI,PK,PM, ■ ■ ■ ,Pm+2)f 
3,1 

xJ^Hpi^ ■ ■ ■ ,PI,PK,PM, ■ ■ ■ ,Pm+2) 

+ X ^klm X%K \Mm{PU ■ ■ ■ ,PI,PK,PM, ■ ■ ■ ,Pm+2)P 
3,1 

x4rnPl> • • • ,PI,PK,PM, ■ ■ ■ ,Pm+2) 



da 



S,d 

NNLO 



-M ^ d^rn+2iPl, ■ ■ ■ ,Pm+2;q)-^ 



m+2 



m+2 



X ^iJk ^'nop \Mm{pi, ...,PI,PK,--- ,PN,PP, ■ ■ ■ ,Pm+2) 



3,0 



J^Hpi^ ■ ■ ■ ,PI,PK, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+2) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



Again, the original momenta of the (m-|-2)-parton phase space are denoted by while 
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Figure 3: Illustration of NNLO antenna factorisation representing the factorisation of both the 
one-loop "squared" matrix elements (represented by the white blob) when the unresolved particles 
are colour connected. 



the combined momenta obtained from a phase space mapping are labelled by /, J, . . .. Only 
the combined momenta appear in the jet function. The phase space mappings appropriate 
to the different cases are described in detail in section 3 below. X^j^i is a four-parton 
antenna function, containing all configurations where partons j and k are unresolved, 
while x^ji^ is a three-parton sub-antenna function containing only limits where parton j 
is unresolved with respect to parton i, but not limits where parton j is unresolved with 
respect to parton k. The factorisation of the phase space is analogous to the factorisation 
at NLO (2.6), such that integration of these antenna functions over the antenna phase 
space amounts to inclusive three-particle or four-particle integrals [22]. 

In single unresolved limits, the one-loop cross section da^]^^Q is described by the sum 
of two terms [44]: a tree-level splitting function times a one-loop cross section and a one- 
loop splitting function times a tree-level cross section. Consequently, the one-loop single 

unresolved subtraction term dcr^j^^^ is constructed from tree-level and one-loop three- 

V S 1 

parton antenna functions, as sketched in Figure 3. Several other terms in dtr^y^^Q cancel 
with the results from the integration of terms in the double real radiation subtraction 
term da^^^Q over the phase space appropriate to one of the unresolved partons, thus 
ensuring the cancellation of all explicit infrared poles in the difference da^^j^Q — dcrj^^^Q. 
Explicitly, the one-loop single unresolved subtraction term is given by the sum of the three 
following contributions: 



X] ~'^ijki^ik) \Mm+l{Pl, ■ ■ ■ ,Pi,Pk, ■ ■ ■ ,Pm+l)f 
X J^+^)(pi, ...,Pi,Pk,... ,Pm+l) 



ik 



dc^iV^LO = m ^^m+l{Pl, ■ ■ ■ ,Pm+i;q)-^ 



(2.12) 



m+l 



m+l 



E 



^iik l-^m(Pl> ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l)|^ J^\P1, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l) 
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+ ^ijk \Mm{Pl, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l)f J^Hpi, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+l) 

(2.13) 

do-So = X] d$m+i(pi, ■ ■ ■ ,Pm+i;q)-^ 



m+1 '^"^+1 



X] '^ijki^ik) X^°op \Mm{Pl, ■ ■ ■ ,Pi,Pk, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+l) 



ik 



X, 



J^\P1, ■ ■ ■ ,Pi,Pk, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+l) 



(2.14) 



In here, X^j^ denotes a one-loop three-parton antenna function. 

Finally, all remaining terms in dafjj^j^Q and do'^ff];^o have to be integrated over the 
four-parton and three-parton antenna phase spaces. After integration, the infrared poles 
are rendered explicit and cancel with the infrared pole terms in the two-loop squared matrix 
element da^'^^^. 

The subtraction terms dcr^^Q, dcr^^^Q and da^^]^Q require three different types of 
antenna functions corresponding to the different pairs of hard partons forming the antenna: 
quark-antiquark, quark-gluon and gluon-gluon antenna functions. In the past [39,40], NLO 
antenna functions were constructed by imposing definite properties in all single unresolved 
limits (two collinear limits and one soft limit for each antenna). This procedure turns out 
to be impractical at NNLO, where each antenna function must have definite behaviours in 
a large number of single and double unresolved limits. Instead, we derived these antenna 
functions in a systematic manner from physical matrix elements known to possess the 
correct limits. The quark-antiquark antenna functions can be obtained directly from the 
e^e~ 2j real radiation corrections at NLO and NNLO [29]. For quark-gluon and gluon- 
gluon antenna functions, effective Lagrangians [45,46] are used to obtain tree-level processes 
yielding a quark-gluon or gluon-gluon final state. The antenna functions are then derived 
from the real radiation corrections to these processes. Quark-gluon antenna functions were 
derived [47] from the purely QCD (i.e. non-supersymmetric) NLO and NNLO corrections 
to the decay of a heavy neutralino into a massless gluino plus partons [45], while gluon- 
gluon antenna functions [48] result from the QCD corrections to Higgs boson decay into 
partons [46]. 

All tree-level three-parton and four-parton antenna functions and three-parton one- 
loop antenna functions are listed in [32]. Their integration over the antenna phase space 
amounts to performing inclusive infrared-divergent three-parton and four-parton phase 
space integrals. Techniques for evaluating these integrals are described in [22,49], and all 
integrated antenna functions are documented in [32]. 

3. Phase space mappings 

The subtraction terms for single and double unresolved configurations described in the 
previous section involve the mapping of the momenta appearing in the antenna functions 
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into combined momenta, which appear in the reduced matrix elements and in the jet 
function. 

At NLO, one needs momentum mappings from three partons to two partons, 

At NNLO, several further mappings are needed. First and foremost, one needs momentum 
mappings from four partons to two partons, 

J^^^-^"^^ : {pi,pj,pk,Pi} {pi,Pl}- 

In addition, repeated mappings from three partons to two partons are also required. 

For the subtraction and phase space factorisation to work correctly, all mappings must 
satisfy the following requirements (specified here for the example of J^^^~*'^^ ) : 

1. momentum conservation: pi + pl = Pi + Pj + Pk + Pi- 

2. the new momenta should be on-shell: pj = 0, p| = 0. 

3. the new momenta should reduce to the appropriate original momenta in the exact 
singular limits, e.g. pi = pi +pj + Pk, PL = Pi in the (i, j, k) triple collinear limit. 

4. the mapping should not introduce spurious singularities. 

The momentum mappings we use follow largely those worked out in [40,43]. The 
different types of mappings needed for our calculation are described in detail in the following 
subsections. 

3.1 Mapping for single unresolved configurations 

In the single unresolved limit where parton j becomes unresolved and i,k are the hard 
radiators, the momenta of the partons k are mapped to pj = {ij) and px = {kj) in the 
following way: 

PI = xpi + rpj + zpk 
PK = {I- x)pi + (1 - r)pj + z)pk , 



where 



' {1 + p)sijk-2rsjk 



z = 



"^{Sij + Sik 

1 

9/, , , X {1- p)Sijk-2rS^J 
z[Sjk + Sik) 



Sijk^ik 



The parameter r can be chosen conveniently, we use [40] 

Sjk 



r = 

Sij + Sjk 
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3.2 Mappings from five partons to three partons 

In the construction of the subtraction terms for the fivc-parton channel, one encounters, 
both the four-parton antenna functions (double unresolved limits), and products of two 
three-parton antenna functions, which compensate for the single unresolved limits of the 
double unresolved subtraction terms. The former require a 4 — >^ 2 mapping, while the latter 
require, in general, a 5 — ^ 3 mapping affecting all momenta. Both types of mappings are 
described in the following. 

3.2.1 Double unresolved configurations 

In the double unresolved limit where partons 12 and become unresolved and ii,i4 are 
the hard radiators, the partons ii, . . . ,Z4 are mapped to the partons ji,j2 with momenta 



Pn = {hhh) , Pj2 = (m«3«2) ■ 



(3.2) 



We will use the shorthand notation (3.2) extensively in the following, keeping in mind 
however that pj-^ and pj^ are linear combinations of all four original momenta with a 
mapping . {pi^,pi^,pi^,pij {Pji,Pj2} given by: 



Pji = xpij^+ ri pi2 + r2 Pis + zpi^ 

Pj2 = (1 - x)Ph + (1 - n)Pi2 + (1 - r2)Pis + (1 - z)Pi4 ■ 
Defining sj^i = {pi^. +Pi,)'^, the coefficients are given by [43] 

S23 + S2A 



(3.3) 



n 

r2 
X = 



\{u, V, w) 



512 + S23 + S24 

513 + S23 + S34 

1 



(1 + p) S1234 



2(S12 + Si3 + S14) 

-ri (s23 + 2 S24) - r2 (s23 + 2 S34) 

+ (ri - r2)- 
1 



Sl4 



Z = 



(1 - p) S1234 



2(si4 + S24 + S34) 

-ri (s23 + 2 S12) - r2 (s23 + 2 S13) 

/ X SI2S34 — ■SI3S24" 

-{ri - r2) 

Sl4 J 

(ri — r2)^ 

1 H o o A(si2 S34, Si4 S23, -513 S24) 



„2 2 
*14 *1234 



+ 



(2 (n (1 - r2) + r2(l - ri))(si2S34 + S13S24 - S23S14) 



su S1234 



+ 4ri (1 - ri)si2S24 + 4r2 (1 - r2) S13S34I] ^ , 
+ + w'^ — 2(uv + uw + vw) . 
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12 unresolved 
*4 ^4 radiator: 

^5 X0(zi,Z2,i4) 




I2 unresolved 
li, radiators 



-32 = hk 
"is = k 



Figure 4: ^g^^^^'^'' i2, «3, *4, «5; ji, j2, J3): both combined partons (^i and ^3) are radiators 
the second step, original parton becomes unresolved. 




12 unresolved 
*4 radiator: 

^5 X0(^l,^2,^3) 




lll2 



■ i3l2 



k unresolved 
k, k radiators 




J3 = k 



Figure 5: (ii, i2, *3, *4, *5; Ji, J2, ^3): one of the combined partons (here /g) becomes un- 



resolved between li and I2 in the second step. 



3.2.2 Iteration of single unresolved configurations 

The four-particle antennae X^-^^ contain by construction all colour-connected double un- 
resolved limits of the (m + 2)-parton matrix element where partons j and k become un- 
resolved. However, X^-^j^ can also become singular in single unresolved limits, where it 
does not coincide with limits of the matrix element. These limits have to be subtracted as 
indicated in eq. 2.9 and as described in section 2.3 of [32]. In the limit where j becomes 
unresolved between i and fe, the four-particle antenna collapses to the three-particle an- 
tenna X^{i,j,k) and a three-particle "remainder matrix element" X^{I, K,l), which also 
has the form of a three-particle antenna. If one of the partons /, K, I becomes unresolved in 
a second single emission, the resulting momenta in these limits coincide with those from a 
double unresolved configuration as defined in subsection 3.2.1. Therefore we need momen- 
tum mappings corresponding to such a "two-step" emission in order to be able to subtract 
the spurious singularities of the four-particle antennae X^jf^^. 

As explained in section 2.2 above and in [32], we have to distinguish between colour- 
connected unresolved partons and almost colour- unconnected unresolved partons. If the 
unresolved partons are colour-connected, we use the momentum mappings ■ 

{Pii ) Pi2 5 Pi:i ' Pit ; Pis } {Ph-,Ph-,Pi3-,Ph} {Pji^Pj2^Pj3}i where the different types B,C 
are described in more detail below. These two-step mappings first map the four partons 
k, . . . ,i4 making up the four-particle antenna to three "intermediate" partons k,k,hj and 
then map the three intermediate partons to two partons ji , j2 ■ The fifth parton 15 = I4 = js 
only acts as a spectator. An additional type of mapping, denoted by J^^~*^~*^\ is needed 
for subtraction terms of two almost colour-unconnected unresolved partons, as defined in 
eq. 2.10 and in [32], and involves redefinitions of all five initial partons. All three mappings 
are depicted in Figures 4-6. 

In the first step, one of the partons {11,12,^3^^4} becomes unresolved. The momentum 
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12 unresolved 
ii, is radiator: 




I2 unresolved 
^3, li radiators 




Figure 6: *2, «3, «4, «5; ji, j2, ja): *3 is the shared hard radiator, 14 — I2 becomes 

unresolved between ^3 and ^4 in the second step. 



is = li is always unaffected by the mapping. The precise definition of the resulting momenta 
h,h: h depends on the mapping. In the first step of the two-step mapping J^^^^^^^ we 
have 



^1 = xiii+ ri 12 + ziii = 1112 

h = h 

I3 = {1- xi) ii + (1 - ri) ^2 + (1 - zi) ii = ui2 

h = h , 



(3.4) 



while for the mappings J^', 



(5~>4^3) 



C 



and JF, 



(5->4-.3) 



K 



, the first step is of the form 



h = xiii + n i2 + zi is = iii2 

12 = 14 

13 = {1- xi) ii + (1 - ri) i2 + (1 - zi) is = 1^12 
h = h ■ 



(3.5) 



The momentum 12 is always the unresolved one, denoted generically by iu in the following. 
The two hard radiators denoted by ia and ib depend on the mapping. In J-l^^^^^\ ii 



and u are the hard radiators, while in ^^^^^^^^ and J^^^'*~^'^' , ii and is are the hard 
radiators. 

Note that mappings of type C apply to all configurations where one of the combined 
partons becomes unresolved in the second step, so they also include the cases where the 
roles of li and 1^ are interchanged. 

The coefficients ri,xi,zi appearing in eqs. (3.4) and (3.5) are given by eq. (3.1), 



,(5^4-^3) 



n 

Xl 
Zl 

p 

2 



+ Sub 
1 

2(Sat( + Sab) 
1 



2(s«b + Sab) L 

4ri(l - ri) SauSub 



(1 + p) Saub -2ri Sub 
(1 - P) Saub - 2ri Sau 



1 + 



SabSaub 



(3.6) 



where now Sub = (i« + ife) etc 
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First step 



Second step 



A/TriDninp" t.vnp 


7 


''ai ''b 


j 


I 


I L 

''ai "'0 


I 


B 




k,k 




h 




u 


C 








h or ^1 


(/i or /s), h 


k 


K 




k,i3 


k 


h 


hi h 


h 



Table 1: Identification of the unresolved, radiator and spectator momenta for both steps of the 
momentum mappings ^^cii*^^ • 

In the second step, one of the intermediate partons {Zi, Z2, 13} becomes unresolved. The 
resulting momenta ji , j2 , js are defined as 



Jl = X2la + r2 lu + Z2h = LL 

h = (1 - X2) la + {l- 7-2) l„ + (I - Z2) k = kk 
h = Ir 1 



(3.7) 
(3.8) 



where again denotes the momentum which is unresolved in the second step, laih ^^re 
the radiators and 1^ does not take part in the second recombination step. The coefficients 
r2,X2,Z2 are defined analogously to eq. (3.6), where now Sub = {}u + h)"^ etc. 

The combination of antenna functions associated with the repeated unresolved singu- 
larity is thus 

X\ia,iuiib)Y\lailuilb) (3.9) 

as indicated in figures 4-6 where and gencrically stand for three-particle tree an- 
tenna functions. The identification of the radiator and unresolved momenta ia, ■ ■ ■ ,lb for 
the different mappings .^b^c^k"^^^ can be read off from Table I and is also illustrated in 
Figures 4-6. 



3.3 Decomposition of antenna functions into sub-antennae 

The antenna phase space mappings require two uniquely identified hard radiator momenta 
and an ordered emission of the unresolved partons. 

With the antenna functions of [32], it is not always possible to uniquely identify the 
hard momenta, especially if more than one final state parton is a gluon. Moreover, in the 
four-parton antenna functions (containing two unresolved partons) at subleading colour, 
the emission is not colour-ordered. It was already outlined in [32] that in both these cases, 
a further decomposition of the antenna functions into different sub-antennae configurations 
is required. 

The tree- level thrcc-parton antenna functions D^{lq,3g,4:g) (quark-gluon-gluon) and 
F^{lg,2g,3g) (gluou-gluon-gluou) contain more than one antenna configuration, since each 
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gluon can become unresolved. Their decomposition was discussed in [32], it reads: 



L>°(1, 3, 4) = 4(1, 3, 4) + 4(1, 4, 3) , (3.10) 
F§{1, 2, 3) = /3°(1, 3, 2) + fl{3, 2, 1) + fl{2, 1, 3) , (3.11) 



with 



4{h 3, 4) = J- f + ^H£5i±ii + fH£H + 5 + Is,,] + Die) ,(3.12) 



'134 



flil, 3, 2) = J- ( 2^ + £12^ + £12^ + 8 \ ^ ^3_^3^ 

«123 V *13S23 S23 ^13 3 j 

With this decomposition, the sub-antennae d^{i,j,k) and f^{i,j,k) contain only a soft 
singularity associated with gluon j, and collinear singularities i \\ j and k \\ j, such that 
i and k can be identified as hard radiators. Soft singularities associated with i on k and 
the collinear singularity i \\ k, which were present in the full antenna functions, are now 
contained in different sub-antennae, obtained by permutations of the momenta. Therefore, 
each sub-antenna can have a unique phase space mapping k) {ij, kj). 

The one-loop three-parton antenna functions Dl{lq,3g,4:g) (quark-gluon-gluon) and 
^3^(1^, 2g, 3g) (gluon-gluon-gluon) can be decomposed according to the same pattern, ex- 
ploiting the fact that each can be written as a function proportional to its tree-level coun- 
terpart plus a function which is not singular in any unresolved limit. 

The decomposition of tree-level four-parton antenna functions is more involved, espe- 
cially since both single and double unresolved limits have to be accounted for properly. It 
turns out to be very useful to introduce the following combinations of three-parton antenna 
functions: 



(1,3,2) = 4(1,3,2) - Al{l,3,2) , 



i?!](l,3,2) = Q0(i^3^2)-Q0(1,2,3) 

Siil, 3, 2) = gO(l, 3, 2) + gO(l, 2, 3) + ^3°(1, 3, 2) . (3.14) 

None of these contains any soft limit or collinear 1 || i limit. Only and contain a 
2 II 3 limit, while is also finite in this limit, owing to the Af = 1 supersymmetry relation 
between the tree-level splitting functions, 

P,g^G{z) + Pgg-.G{Z) = Pgg-.Q{z) + Pgg^Q{l - z) . (3.15) 

Among the tree-level four-parton antenna functions, only Al{lq,3g,4:g,2q) (quark- 
gluon-gluon- antiquark at subleading colour), 1^4(15, 3^, 4^, 5^) (quar k-gluon-gluon- gluon) , 
E^{lq, 3qi,Aqi,bg) (quark-quark-antiquark-gluon at leading colour), F4 (1^, 2g, 3g, 4g) (gluon- 
gluon-gluon-gluon) as well as G\{lg, 3g, 4g, 2^) and G\{lg,3q, 4g, 2g) (gluon-quark-antiquark- 
gluon at leading and subleading colour) must be decomposed into sub-antennae. In the 
context of the three-jet calculation discussed here, F^, G\ and G\ do not contribute and 
will not be discussed further. 



-16- 



The decomposition of ^4(1^,3^,4^,2^) is needed because both gluons can become 
colhnear either with quark Ig or with antiquark 2q. The two possible phase space mappings 
are of the ^(5"*^) type described in section 3.2.1 

(a): (1, 3, 4, 2) ^ (134, 243) (b): (1, 4, 3, 2) ^ (143, 234) 

each allow only gluons adjacent to quark or antiquark to become collinear. To disentangle 
the different sub-antennae, it is sufficient to partial-fraction the antenna function in the 
different collinear denominators, as done in [32]. With these, we decompose 

Al{l, 3, 4, 2) = i° „(1, 3, 4, 2) + A%{1, 3, 4, 2), (3.16) 

with 

io,(l,3,4,2) = ^0(1, 3, 4, 2) +0^(2,4,3,1) , 

10^(1,3,4,2) = i" , (1,4, 3, 2) , (3.17) 

where a'l{i,j, k, I) contains only singularities for i \\ j or A; || / and was defined in [32]. With 
this decomposition, ^4^(1,3,4,2) contains only 1 || 3 and 2 || 4 singularities, and can be 
used with the (1, 3, 4, 2) (134, 243) mapping. Since this decomposition is straightforward, 
we refrain from spelling it out explicitly in the subtraction terms presented in subsequent 
sections. 

The leading-colour quark-quark-antiquark-gluon antenna function £^2(1^, 3g', 4^/, 5^) 
contains limits where cither the quark-antiquark pair (3g',4g/) or the gluon 5^ can become 
soft. Since these limits yield different hard radiator partons, they can not be accounted for 
in a single phase space mapping, but require two separate ^(^^^) mappings: 

(a): (1, 3, 4, 5) ^ (134, 543) , (b): (1, 5, 4, 3) ^ (154, 345) . 

By analysing the different triple and single collinear limits of E^, one finds the following 
decomposition: 

El^{l,3,4,5) = b2{1,3,4,5) + El{5,4,^)Ql{hm,m) , (3-18) 
El,{l, 3, 4, 5) = e2{1, 3, 4, 5) - E°„(l, 3, 4, 5) . (3.19) 

After this decomposition, E^^^ can be used with mapping (a) and E^ ^ with mapping (b) . 
The above decomposition also ensures a well-defined behaviour in all double and single 
unresolved limits (see [32] for a definition of the splitting factors): 

£;4%(1,3,4,5) ' 5i5(3,4), 

„ 3„/||4^,,5„^0 1 

El,{l, 3, 4, 5) " SvM^) — P,^^g{z) , 

i^4%(l,3,4,5) ''"M''-' P^^l-^^-'iw,x,y), 
£;2,a(l,3,4,5)+£;2,6(l,3,4,5) """'V' Pus^G{w,x,y) , 

£;°,(1,3,4,5) ^"^"^ p^^^q(^) p^-_^^(y) , (3.20) 
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3«/|IV 1 



S34 



P,,-_^GWd^(l,(34),5)+ang. , 



El,{l,3,A,5) — Pqq^Giz) d^(l,5,(34)) +ang. , 

£;4%(1,3,4,5) '-^^ —P,,^Qiz)El{{15),3,4), 

E%{1,3A,5) —Pgg^Q{z)El{l,3,{45)), 
S45 



(3.21) 



while all other limits are zero. It can be seen that only the triple collinear 3 || 4 || 5 and the 
single collinear 3 || 4 limits receive contributions from both phase space mappings. This is 
unavoidable, since both these limits match onto the double soft (3g/,4^/) and the soft 
5g limits, which belong to different mappings. 

The decomposition of the quark-gluon-gluon-gluon antenna function 0^(1^, 3^,4^,5^), 
is more involved, since any pair of two gluons can become soft. We consider four different 
jr{5-*3) niappings: 



(a): (1,3,4,5) ^(134,543), 
(c): (1,3,5,4) ^(135,453), 



(b): (1,5,4,3) ^(154,345), 
(d): (1,5,3,4) ^ (153,435). 



The numerous different double and single unresolved limits of this antenna function can be 
disentangled very elegantly by repeatedly exploiting the M = 1 supcrsymmctry relation [37] 
among the different triple collinear splitting functions [35,37,50]. Using this relation, one 
can show that the following left-over combination is finite in all single unresolved and 
double unresolved limits: 



L>0/(1,3,4,5) = ^2(1, 3, 4, 5) 



A0(1,3,4,5)+A0(1,5,4,3) 



~ (^2(1,3,5,4) +£:0(1, 5, 3, 4)) +i°(l,3,5,4) 

-£^2(1, 5, 4, 3) + S2(l, 5, 4, 3) + C2{1, 4, 5, 3) 

-E2{1, 3, 4, 5) +£2(M, 4, 5) + C4°(l, 4, 3, 5) 

+^0(1, 3, 4) Sl{{13), (43), 5) + A^l, 5, 4) Sl{{15), (45), 3) 

+^0(3, 4,5)50(1,(54), (34)) 



(3.22) 



Starting from the terms in this expression, the following sub-antennae can be constructed: 

i^o„(l,3,4,5) = ^i^0/l,3,4,5)+A0(l,3,4,5)-^40(l,3,5,4) 
+^0(l,3,4)530(('l3),(43),5) 

+1^0(3,4, 5) [Slil, (M), ('34)) - i?0(l, (54), (U))) 
-El{5, 4, 3)Q°(1, (34), (M)) - yl°(l, 3, 4) eIOU), (43), 5) 
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-A0(1,3,4)Q0((13),5, (43)), 
Z)°fe(l,3,4,5) = A%(1,5,4,3) , 

Dl,{l, 3, 4, 5) = 3, 5, 4) - e2{1, 5, 4, 3) + 5, 4, 3) + C^ih 4, 5, 3) 

+^0(3, 4, 5)Q!1(M54), (34)) + 3, 4) i^^O^^ 

+aO (1, 3, 4) ((13) , 5, (43) ) + (4, 5, 1) ((15) , 3, (45) ) , 

L>o^(l,3,4,5) =i^yi,5,4,3). (3.23) 

With this decomposition, each D^j^ contains only singularities appropriate to phase space 
mapping (i). The sum of the adds to D^: 

Dl, + D% + Dl, + Dl, = Dl , (3.24) 

such that only must be integrated analytically over the antenna phase space. 

The above decomposition disentangles the different double and single unresolved limits: 



i^4,a(l,3,4,5 



D 



4,& 



(1,3,4,5)+!)°^ 
D0,(1,3,4,5)+D4% 

L>4%(1, 3, 4, 5) + i^4%(l, 3, 4, 5) + 1^4% 



D 



4,6 



D4%(l,3,4,5)+i?4% 



^0 (1,3,4,5)+Z)0 (1,3,4,5)+Z)0 



-^4,& 



i^4° (1,3,4,5) +1^4% 



Z)0„(l,3,4,5)+i?4% 
Z^o,(l,3,4,5)+Z)°rf 



1,3,4,5 

1,3,4,5 
1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 
1,3,4,5 
1,3,4,5 
1,3,4,5 
1,3,4,5 

1,3,4,5 



4g- >0,5g ->0 

3g^0,5g->0 

l^ll Sg.Sg ^O 
4sl|59,3g->0 



S'l345 , 
5'l543 , 



5'l34 'S'l54 , 

'S'4;315(^) — Pqg^qi} " z) , 

'S'l;345(2;) Pgg^aiz) , 

S45 



1,1139,49^0 1 

^ 'S5;43l(2) Pqg^Q{z) 

Sl3 



lq|| 5g,4g -^0 

1,1139,59^0 



'S'3;45l(2:) Pqg^qiz) , 

Sl5 

'S'4;513(^) Pqg^qi'^ " ^) , 

Sl3 



391149,59-^0 / N 1 D / N 

>• '5i;543(2;) ^39^0(2:) , 

S34 

lgl|3g||49 
I9 II59II 49 
I9 II39II 59 

— ' P345^G[W,X,y) 

1,1139,491159 1 



^i34-^Q(^«,a;,y) , 
Pi54-^Q{'w,x,y) , 

Pi35^Q{w,x,y) , 



S13S45 

S15S34 



Pqg^qi^) Pgg^aiy) , 

Pqg^qiz) Pgg^ciy) , 

(3.25) 



D2„(l,3,4,5)'^°5i34d^(l,4,5), 
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L>°c(l,3,4,5)+L>°^ 



-^4,0 

-'^4,6 
-^^4,6 



i54%(l,3,4,5)+i^4% 



-^4,a 



4,c 



-^4,6 



-^4,(i 



-^4,a 



i?4%(l,3,4,5)+Z)4% 



4,6 



^4%(1,3,4,5)+D0 



1,3,4,5 
1,3,4,5 
1,3,4,5 
1,3,4,5 
1,3,4,5 
1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 

1,3,4,5 



5i34c/^(l,5,4) 

'-^^ 5345 ^^(1,3,5) 

^11!° 5345^0(1,5,3) 

5154^0(1,4,3) 

5x54^0(1,3,4) 
-^-P,,_Q(z)dO((13),4,5), 

Sl3 



1 

Sl3 



P,,^Q(z)dO((13),5,4) , 



— P,,_Q(z)dO((15),4,3), 
-^P,,^Q(z)dO((15),3,4), 

— Pgg^Giz) dg(l, (34), 5) + ang. , 

•S34 

— Pgg^Giz) 4ih 5, (34)) + ang. , 
S34 

59 1 n /.A jO 



S45 

5a 1 



P3,^G(2;)c/^(l,(45),3)+ang. , 



? — Pgg^aiz) c?^(l, 3, (45)) + ang. . (3.26) 
S45 

All other limits are vanishing. It can be seen that certain limits are shared among several 
antenna functions, which can be largely understood due to two reasons: 

1. in a gluon-ghion collincar splitting, either ghion can become soft, and the gluon-gluon 
splitting function is always shared between two sub-antennae, as in (3.10), (3.11) to 
disentangle the two soft limits. 

2. the unresolved emission of gluons 3^ and 5^ is shared between the mappings (c) and 
(d) according to the decomposition of the non-ordered antenna function ^4, which 
distributes the soft limit of either gluon between both mappings. 

3.4 Angular terms 

The angular terms in the single unresolved limits are associated with a gluon splitting into 
two gluons or into a quark-antiquark pair. They average to zero after integration over the 
antenna phase space. To ensure numerical stability and reliability, this average has to take 
place within each phase space mapping. We have checked this to be the case for the above 
decompositions of and D\. The angular average in single collinear limits can be made 
using the standard momentum parametrisation [10,51] for the i \\ j limit: 



with 2pi ■ pj 



z 2p - n ^ 



Pi 



(1 - z)pf' - fe^ 



1 — z2p ■ n ' 



z{l - z) 



p2 = ^2 ^ Q 
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In this denotes the cohinear momentum direction, and is an auxiliary vector. The 
collinear limit is approached by A;^ — > 0. 

In the simple collinear i \\ j limit of the four-parton antenna functions X2{i,j,k,l), 
one chooses n = pkto be one of the non-collinear momenta, such that the antenna function 
can be expressed in terms of p, n, k± and pi. Expanding in A;^ yields only non- vanishing 
scalar products of the form pi ■ k±. Expressing the integral over the antenna phase space 
in the {p, n) centre-of-mass frame, the angular average can be carried out as 



^- r d</> {p, ■k^)=o, ^ r dct> {pi . k^f = -ki piPiniiPi . (3.27) 

27r Jq 27r Jo P-n 



Higher powers of k'^ are not sufficiently singular to contribute to the collinear limit. Using 
the above average, we could analytically verify the cancellation of angular terms within 
each single phase space mapping, which is independent on the choice of the reference vector 

The remainder of this section has been modified compared to the original version of the 
paper: In the N'^ and N'^ colour factor, the angular averaging is not sufficient to cancel 
the 1/e poles in the four-parton one loop subtraction terms [B]. In either of these colour 
factors, the difference dcrjj^^^Q — da^^]^Q contains left-over poles of the form 

1x0(1, i, 2) Y,\U,j, 2l) jf (fz, j, 2i) [sz^ + sg. - s^J - s^J - s^^, + } , (3.28) 

where X3 and are tree-level three-parton antenna functions. Contrary to statements 
made in [B] , these terms do not appear in the colour factor Np N in our implementation. 

Furthermore, for these two colour factors the five-parton subtraction terms themselves 
do introduce spurious limits from large angle soft radiation. The single soft limit of i or /c 
in (6.3) is non- vanishing. Instead, it yields (soft i): 



+ ldl{2,k,j)Al{l,jk,2k) [S^~^ + S^^^^-S^^^^^-Suj-Su2 + S,,,]+{1^2) 
- 1^0(1, k, 2)Alm),j, (2k)) [5,^. . + 5,^. . - 5,^, - Su, - S,,, + (3.29) 



2 

with 



Sahc = (3.30) 

To account for this large angle soft radiation, a new subtraction term dcr^^y^Q is intro- 
duced. This term is added to the five-parton subtraction term dcr^^^^Q, and its integrated 

V S 1 

form is subtracted from the four-parton subtraction term dd^y^^Q, cancelling the left - 
over 1/e terms and adding new finite contributions to the four-parton and the five-parton 
subtraction term. 

The new subtraction term du^^^^Q contributes only in the N"^ and N'^ colour factors. 
Its contribution to N"^ reads: 

^'^iNLO,N^ = ATg Ar2 d$5 (pi , . . . , P5 ; g) 1 ^ J 

■ (i,i,fc)GPc(3,4,5) I 
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+ 



^0 



' Q Q -I- Q Q 

'{{li)k)i{{ji)k) 2i{{ji)k) 2iUi) 2j((H)fe) 

(3), 



4-—( Q — Q — . <? _(_ Q . C I o A 

2V ilk)k{jk) 2k{{jk)i) 2k{jk) 2k{{lk)i) 2k{\k)) 



X 



ig, {jk)g) Almk)i)^, {{jk)i)g, 2g) {p^i,^i,P2) 





^2V ((2i)fc)»((r0fc) 



Sr. 



C _|_ Q ~ _ C ^ i_ C 

li{{ji)k) ^ li((2i)fc) ^ H(2j) 



(2i)j(ii) 

xd°((2r)^-, /c,, c^^) ^0(1,, (c^)^, (^)^-) jf\pi,^k,m)k) 

I I/' c.^^ _ . c I q q |_ q 

2\ {{2k)i)k{{jk)i) {2k)k{jk) lk{{jk)i)^ lk{jk) lfc((2fc)i) " 

x4{{2k)^,ig, (jk)^) Al(lg, ((jk)i)g, ((21^)-) jf\vi,nM'.^)i. 

~ 2 ('^({n)fe')i((2i)fe) ~ " ^{^2i)k)ij ^(M)V ^ ^W)i3 ~ ^Jli)i{2i) 

x^O((H)^, kg, (2i)^-) Al{((li)k)^,jg, ((2i)k)g) 4'\^k,Pj,m)k) 



l(s 

2V ((lfc)i)fc((2fc)i) 



^((lk)i)kj ^((2k)i)kj ^ ^{lk)kj + ^^(2^)^^ ^{lk)k{2k) 

(3), ^ 



xAl{{lk)^,ig,{2k)^)Almk)i)^,jg,{{2k)i)^)4'\^^^^ (3.31) 



The new contribution to the five-parton subtraction term is: 



\ E ( 

(i,j)6(3,4) 



S——- -\- s —S' 

((W)5);j {{2i)b)ij ((H)5)i((2i)5) 



x^O((li)^,5,,(2i)^-)^0(((H)5)^,j,,((2i)5)^0#(j^5,]^5,P^^ 



o) 



(3.32) 



These large-angle soft subtraction terms contain soft antenna functions of the form Sajc 
which is simply the eikonal factor for a soft gluon j emitted between hard partons a and c. 
Those soft factors are associated with an antenna phase space mapping {i,j,k) {I,K). 
The hard momenta a, c do not need to be equal to the hard momenta i, k in the antenna 
phase space - they can be arbitrary on-shell momenta. 

The integral of each of these soft antenna functions over the antenna phase space can 
be written as 



c;ik — J ^^Xijk 



s. 



ajc 



r^(l - e)e^^ 
r(l - 3e) 



where we have defined 



Xac,IK 



h In [xaciK) + e Li2 

e V XacJK 



SacSlK 



{Sal + SaK){ScI + Sck) 



(3.33) 
(3.34) 
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LO 


7* 


^ qqg 


tree level 


NLO 


1 


— > a no 


one loop 




7* 




tree level 




7* 




tree level 


NNLO 


7* 


^ qqg 


two loop 




7* 


qqgg 


one loop 




7* 


QQQQ 


one loop 




7 


qqqqg 


tree level 




7 


qqggg 


tree level 




7 


ggg 


(one loop)^ 



Table 2: The partonic channels contributing to e+e 3 jets. 
So that the integration of the new N"^ subtraction terms reads 

J ^'^x,j,da^NLO,N^=^^^^ (^) d$4(pi,...,P4;g) X 

(ij)e(3,4) l 

(j5,,2,-)jf (pri,p7i,P2) + (1 - 2) 

yl°(l,,i,,2,-)^0((li)^,j,,(2i)-)jf (3.35) 

while for the A'^'^ term the integration of the 5-parton contribution over the antenna phase 
space yields 

j d$x,,,d(7^jv^o^^o =iV4ArO d^4{pi,P2,P3,P4;q) 

\ ('^(Ii)j;12 + '^(5^)j;12 - '^(II)(5^;12 - '^li;12 - 52J;12 + 5^^^ (3.36) 

(jj)e(3,4) 

4. Pcirton-level contributions to e+e~ — > 3 jets up to NNLO 

Three-jet production at tree-level is induced by the decay of a virtual photon (or other 
neutral gauge boson) into a quark-antiquark-gluon final state. At higher orders, this process 
receives corrections from extra real or virtual particles. The individual partonic channels 
that contribute through to NNLO arc shown in Tabic 4. 

According to the structure of the coupling to the external vector boson, one distin- 
guishes non-singlet and singlet contributions. The non-singlet contributions arise from 
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the interference of amplitudes where the external gauge boson couples to the same quark 
lines, while the pure singlet contribution is due to the interference of amplitudes where 
the external gauge boson couples to different quark lines. Up to NLO, only non-singlet 
contributions appear. It is only at NNLO, that the first non-vanishing singlet terms are 
allowed. These appear in the tree-level 7* qqqqg process, the one-loop 7* qqgg 
and 7* qqqq processes and the two-loop 'j* ^ qqg process. All these processes yield 
both non-singlet and singlet contributions. The 7* ggg process, which is mediated by a 
closed quark loop, is entirely a singlet contribution. In four-jet observables at ©(a^), the 
singlet contributions were found to be extremely small [52]. Also, the singlet contribution 
from three-gluon final states to three-jet observables was found to be negligible [53]. 

Matrix elements and subtraction terms at NLO and NNLO can be naturally decom- 
posed according to their colour structure. The cross section at NLO receives contributions 
from three different colour factors: 

^CTNLO = daNLO,N + d(TjvLO,l/Ar + ^(^nlo,Nf ■ (4-1) 

The NNLO contribution to the cross section receives contributions from seven different 
colour factors: 

dcTiVATLO = d<^NNLO,N'2 + d<^NNLO,N0 + d(JNNLO,l/N'^ 

+daNNLO,NF N + d(^NNLO,NF/N + d(T NNLO,Nj. + d(^NNLO,NF^^ ■ (4.2) 

The first six terms in this equation are non-singlet contributions, the last term is the 
numerically unimportant singlet contribution. 

In the following, we list the matrix elements for the contributing partonic channels 
shown in Table 4 and discuss their structure. 

4.1 Tree- level matrix elements for up to five partons 

The tree-level amplitude for a virtual photon to produce a quark-antiquark pair 

and (n — 2)-gluons, 

i*{q) q{pi)q{p2)g{P3) ■ ■ ■ g{pn) 

can be expressed as sum over the permutations of the colour ordered amplitude ^ of 
the possible orderings for the gluon colour indices 

K{n-2)g = i<V2gr-' E (^"' • • • ^"'')ni2 •^A,n(Pl,P3, ■ ■ ■ ,Pn,P2) • (4.3) 

(i,...,fe)GP(3,...,n) 

The squared matrix elements for n = 3, . . . , 5, summed over gluon polarisations, but 
excluding symmetry factors for identical particles, are given by, 

\M^g/ = NsA',{lg,3g,2g), (4.4) 



E NAl{l„ig,jg,2g) - liO(l„ 3^,4^,2,-) 



N 

(i,i)GP(3,4) 



(4.5) 



-24- 



where, 



(j,j,fe)eP(3,4,5) 



+ 



iV2 + 1 
iV2 



iV„ = 4;ra5;e^(/)("-^) (iV^ - l) |>1 



^bi^qj 4g, 5g, 2qy 



'99 1 : 



(4.6) 



(4.7) 



and 



K/ = 4(l-6)g^ 



(4.8) 



The squared colour-ordered matrix elements ^43, and are given in [32]. For the five 
parton case [54], 



^5(1?! ^giis! ^ff) 2g) 



qq\ 



^5(195 ^95 Jgj ^95 2q) |A1 



|2 



M%r^{pi,Pi,Pj,Pk,P2) + M\5{pi,Pi,Pk,Pj,P2) + M%ri{pi,Pk,Pi,Pj,P2] 



^bi^Q'h'jg' ^5) 



|2 



X] MA,Q{pi,Pi,Pj,Pk,P2] 
ii,j,k)€P{3,...,5) 



(4.9) 

) 

, (4.10) 
(4.11) 



In the subleading colour contribution ^5, gluon k is effectively photon- like, while in the sub- 
subleading colour contribution (also called Abelian contribution), A^, all three gluons are 
effectively photon-like. Photon-like gluons do not couple to three- and four-gluon vertices, 
and there are no simple collinear limits as any two photon-like gluons become collinear. As 
a consequence, the only colour-connected pair in ^5 are the quark and antiquark. 
The tree-level amplitude for 



1*{q) Q{Pi)Q{P2)q'{P3)q'{P'i) 



is given by 



^qqq'q' = if^lQ^ ^qiq2^q3qi {^111^312 " ;^'^ili2^i3i4^ M%^i{pi,P2,P-i,PA) 

+(1^3,2^4), (4.12) 

where 5q^q,^5q.^q^ indicates the quark flavours. The amplitude M% 4^{pi,P2,P3,P'i) thus de- 
notes the contribution from the gig2-pair coupling to the vector boson. The identical quark 
amplitude is obtained 

(4.13) 



<-..- = <aVa-'-(2--4). 



qqqq ^'^qqq'q- 
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The resulting four-quark squared matrix elements, summed over final state quark 
flavours and including symmetry factors are given by 



^4/ = El^^^wl'+Ei^' 



NpB^ilg, 3g, 2g) - ^ {C^^{lg, 3„ 4g , 2g) + C°(2,-, 4g, 3„ 1,)) 



(4.14) 



where 



C2{lg,3g,4g,2g) I' = -Re(M%Api,P2,P3,P4)M''j^\(.Pi,P4,P3,P2)] , (4.15) 



B2{lg,3g,4g,2g) | = Re ^Al ^,4 (^1 , ^2 , , P4 ) ^ ^4 ^3 , , Pi , ^2 ) Y (4.16) 



Explicit expressions for and C4 are given in [32]. The last term, B^, is proportional to 
the charge weighted sum of the quark flavours, Np^^, which for electromagnetic interactions 
is given by. 



(4.17) 



It is relevant only for observables where the flnal state quark charge can be determined. 
There are four colour structures in the tree-level amplitude for 



l*{<l) q{Pi)<l{P2)q'{P3)q'{Pi)g{Pb 



which reads 



Mqqq'q'g = iSig^ V2Sg^gJg,,g^ 



,P2,P3,P4,P5) - ,P2,P3,P4,P5) 

+^3i2 Ku-M%5(P'L,P2,P3,P4,P5) - ^ T^^'X 1 i 2 S J > P2 , P3 , , ) 

+ (1^3,2^4). (4.18) 

The amplitude ■M^^^^{pi,P2,P3,P4:,P5) for x = a, . . . , d denotes the contribution from the 
Q1Q2 pair coupling to the vector boson. Due to the colour decomposition, the following 
relation holds between the leading and subleading colour amplitudes: 

■A^B^5(Pl,P2,P3,P4,P5) = M^j^'^^{pi,P2,P3,P4,P5) + M'^J',^{pi,P2,P3,P4,P5) 

= ■M.%%iPl^P2,P3,P4,Pb) + ■M!]f^^{pi,P2,P3,P4,Pb) ■ (4.19) 
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As before, the identical quark matrix element is obtained by permuting the antiquark 
momenta, 

= - (2 4). (4.20) 

The squared matrix element, summed over flavours and including symmetry factors is 
given by, 



NNf 5g, 3,., 2,-) + 3,., 5^, 2,-)) 

+^ (Sf (1„ 5„ 2,-; 3,,, %0 + Bl'\l,, 2,-; 3,,, 5„ - 25°'^(1„ 2,-; 3,,, 5,)) 

"^"5(1^, 3g, 4g, 5g, 2g) + ^ ^ ^(75(15, 3g,4g, 5g , 2g) + (75(2^,45, 3g, 5g , Ig)^ 

-iVATf,^ (4°'«(l„55,4,-,;3g,,2g-) +Sf (lg,4,-,;3g,,5g,2,-) - 4,-,; 3,,, 2,-, 5^)) 

H — ^ (-^5'''(l5'5s'25;3g/,4g/) + ^^''^(Ig, 2^; 3^/, 5^, 4^/) + ^^'''(Ig, 2^; 3g/, 4g/; 5^)^ , 

(4.21) 

where for x = a, . . . , e 



^5''(- ■ •) = \M%{pi,P2,Pz.PA.P,)\\ (4.22) 

■ ■) = Re (M^f5(i'l>^'2,P3,P4,P5)A4^f/(p3,P4,Pl,i>2,P5)) , (4.23) 

and 

C0(l„3g,4g-,55,2g-) |mO/ = -2Re(^A4^;'5(pi,p2,i?3,P4,P5)A^^'^(pi,P4,i?3,P2,P5) 

+ si's (^'l ' P2 , P3 , P4 , P5 5^5^ (Pl , P4 , P3 , P2 , P5 ) 
+ -^B"5(f'l' ^^2, P3, P4, P5)-Mb1^^ {P3,P2,Pl , P4, Ps) 

+ A4^|'5(P1,P2,^»3,P4,P5)A4b^5^(P3,P2,^»1,P4,P5)^ , 

(4.24) 

C5{lq,Sg,Ag,5g,2g) | M | ^ = " Rc ^ (^1 , ^2 , ^3 , ^4 , ^5 )-M £^5^ (pi , ^4 , 153 , ]32 , j^S ) ^ • 

(4.25) 

4.2 One-loop matrix elements for up to four partons 

The renormalised one-loop amplitude M^^g for a virtual photon to produce a quark- 
antiquark pair together with a single gluon, 

7* (9) Q{pi)q{P2)g{P3) 
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contains a single colour structure such that 

-2 



Ml,, = ieV2g ( J-^) T^X^X3iPi,P3,P2) 



(4.26) 



Unless stated otherwise, the renormalisation scale is set to = q^. 

The interference of the one-loop amplitude with the three-parton tree-level amplitude 
(4.3) is given by 



2Re (MgM^J = iV3 (g) A^'Hu, 3„ 2,) , 



(4.27) 



where 



4'''°^(lg,3^;,25) = (^N[Al{lq,3g,2g)+Al{si23)Alilg,3g,2q)] 

i^(l„ 3^, 2,-) + ^1(5123)4(15,3^,2,-)] +NFAl{lg,3g,2g)j , 

(4.28) 



1 

'iV 



where A^, A\ and A\ are given up to 0{e^) in [32]. 

Moreover, the one-loop process 7* — > ggg also yields three-jet final states. Since this 
process has no tree-level counterpart, it does only contribute at NNLO. Its amplitude can 
be denoted as [53] 

K99 = (if^) d"'"'"'Mh^^(j>i,P2,Pz) . (4.29) 

The one-loop corrections to 7* — ^ 4 partons have been available for some time [21]. 
The one-loop amplitude for 



7* (9) q{pi)q(p2)g{p3)g(p4.) 



contains two colour structures, 



(i,j)eP(3,4) ^ 

+NfM]i'^^ {Pl ,Pi,Pj,P2)-\ ^ ^^4 (Pi ^PhPi^P2)\ 

+ ^'5"'"^'^ni2A4^J(Pi,^'3,P4,P2) , (4.30) 



where 



MY^^{pi,P^,Pa,P2) = MY^^{pi,pi,p^,p2) 



(4.31) 
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The "squared" matrix element is the interference between the tree-level and one-loop 
amplitudes, 



1199 <1<199\ 



+NNF,^AY{lq,i9,j9,'^l)) 



(i,j)6P(3,4) 



A^ 3y , 4y , 2g) A^ (Ig, 3g, 4g, 2q) -^^A^ (Ig, 3(,, 4^, 2q 



■-^^4 (Igj 3g, 4g, 2q) H jf'^^ (Ig; 4g, 2g) 



(4.32) 



where for x = a, . . . , c?, 

AY{lq,ig,jg,2q) \M%\^ = Re (M\%pi,Pi,Pj,P2)M\]^{pi,Pi,Pj,P2)) , (4.33) 
i^'^ (l„3<„4g,2,-) \M%\^ = Re {mXM^P3,P^,P2)M%{pi,P3,P4,P2)) , (4.34) 

and 

M]i^4^{pi,P3,P4,P2) = M](^^{pi,P3,P4,P2) + M]i^^{pi, P4, Ps, P2) ■ (4.35) 

The renormalised singularity structure of the various contributions can be easily writ- 
ten in terms of the tree-level squared matrix elements multiplied by combinations of infrared 
singularity operators [55], for which we use the notation defined in [32]. Explicitly, we find 

roles{AY{lq,ig,jg,2g)) = 2 (lW(e,SH) + lg(e,Si,) +I^y(e,s,-2)) ^^(l,, i^, j^, 2,-), 

(4.36) 

Voles{AY{l„igJg,2g)) = 24^(6, .si2)A^(lg, ig, jg, 2,~), (4.37) 
Voles{AY{lg,ig,jg,2g)) = 2 (O(6,si0 (e,Si,) +l(^;p(e,s,-2)) ^2(l„i„ j„ 2,-), 

(4.38) 

P0lesiAYilg,Sg,ig,2g)) = 

(2lW(e,S34) +lg(e,si4) +4^(e,S23) + lS';(e,si3) + 4^^ i5(l„ 3„ 4„ 2,-), 

(4.39) 

VolesiAYilg,SgAg,2g)) = 21^^ (c, Si2)i2(l„ 3^, 4^, 2,), (4.40) 

V0les{AYilq,SgAg,'2q)) = 

(21^^(6, .34) + ^14) + l^glA^, S23) + I^fC^, Sn) + l^gl^i^, S2a)) 

xi°(l„35,4^,2,-), (4.41) 

roles{AYilg,^gAg,'2g)) = 

2I^^jie, si2) + 2I«(6, S34) - I«(e, su) - l£^(6, ^23) - lg(e, ^13) - 4^(e, ^24)) 

xAl{lg,3g,Ag,2g). (4.42) 
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As at tree-level, the one-loop amplitude for 

q{pi)Q{p2)q'{P3)q'{p4) 

contains two colour structures, 

^2 



[NMBy{pi,P2-,PZ-,PA) - JzMB^4kPuP2,P^,Pi) + NFMB,i{puP2,PZ,PA) 



+NFM^^^i {pi ,P2,P3,P4) 



+ (1^3,2^4) , 



(4.43) 



where 

M]f^ipi,P2,P3,PA) + MsAipi,P2,P3,P4:) = Mb4{pi,P2,P3, ,P2,P3,P4)- 

(4.44) 

As before, the identical quark matrix element is obtained by permuting the antiquark 
momenta. 



M, 



Miw-(2-4). 



(4.45) 



9959 ^"qqq'q' 

Summing over flavours and including symmetry factors, we find that the "squared" 
matrix element, is given by 

9,9' 



9999 



"^"4' {lq,3q,4q,2q) + j^C[^ (Ig, 3g, 4g, 2g) -^C4''^(lg,3g,4q,2q,) 
— C^''^{2q, 4q, 3q, Iq) + j;^C^''^{2q, 4q, 3q, Iq) —C^'-^{2q, Aq, 3q, Iq) 

+NNF,^Bl'''{lq, 3q,4q, 2q) - ^Bl'\lq, 3qAq, 2q) + iV^iVi.,^^|''^(lg, 3q, Ag, 2q) 



NNpE^ilq, 3q, 4q, 2q) - ^B\\lq, 3„ 4,-, 2,-) + Ar|i?]'^(l„ 3„ 4q, 2q) 



N 



(4.46) 



where for x = a,b,c 



Bl'''{lq,3q,4q,2q)\M 



|2 



s]'"(l„3„4,-,2,-)|7W0/ 



Re [M]^'l^ipi,p2,P3,P4)M'j^l^{pi,p2,P3,P4)j , (4.47) 
R-e (M^BAiPi^P2,P3,P4)M^B4iP3,P4,Pi,P2)) , (4.48) 



and for x = d,e, f 



Cl'^(l„3„4,-,2g-) |A^o/ = -Re (m]^%Pi,P2,P3,P4)MI]^{pi,P4,P3:P2)) • (4.49) 
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Using the infrared singularity operators of [55], we can extract the singular contribu- 
tions of the renormalised one-loop contribution as, 

PoZes(S]'«(l„3„%,2,-)) = 2 (ij^ (e, sm) + ij^ (e, S23)) 3„ 2,-), (4.50) 

rOlesiBl'\lg,3g,ig,2g)) = 

2 (24^ (e, 514) - 24^(6, S13) + 24^(6, S23) - 24^(6, S24) + 4? (e, 512) + 4$ (e, ^34)) 

xS0(l„3„4,-,2,-), (4.51) 
Voles{Cl^''{l„3„4q,2g)) = 2 (4^ (e, ^13) + 4^ (^^ ^24)) C4°(l„ 3^, 4,-, 2,-), (4.52) 

2 (4^ ^12) + I^^^ (e, S14) + I^^^ (e, S23) + I^^^ (e, S34) - 4^ ^is) - 4^^ ^24)) 

xC40(l„3g,4,-,2,-), (4.53) 
Voles{Bl'''{lg,3g,4g,2g)) = 2 (4^)(e,si4) +4? (e,S23)) ^^(l,, 3„ 4,-, 2,-), (4.54) 

PolesiBl'\lg,3q,4g,2g)) = 

2 (24^ (e, su) - 24^^ (e, ^13) + 24^^ (e, S23) - 24^^ (e, ^24) + 4? (^^ ^12) + 4¥ ^34)) 

x540(l„35,4,-,2,-). (4.55) 

4.3 Two-loop matrix elements for three partons 

The renormalised two-loop amplitude M^^g for a virtual photon to produce a quark- 
antiquark pair together with a single gluon, 

i*{q) q{pi)q{p2)9{P3) 

contains a single colour structure such that 

/ 2 \ 2 

M'ggg = ieV2g{^-^j T^%MX3iPi,P3,P2) . (4.56) 

At NNLO, there are two contributions. One from the interference of the two-loop and 
tree-level amplitudes (4.3), the other from the square of the one-loop amplitudes given 
in (4.26). These terms were computed in [18] by reducing the large number of two- loop 
Feynman integrals to a small number of master integrals, using integration- by-parts [56] 
and Lorentz-invariancc [57] identities, solved using the Laporta algorithm [58]. The relevant 
master integrals (two-loop four-point functions with one off-shell leg) were then derived [59] 
from their differential equations [57,60,61]. 

The resulting virtual three-parton contributions are given by, 

2Re (m^M^^J = iV3 (g)' 4'^°)(1„3„2,) , (4.57) 
Re (m^^MiJ = iV3 (^)' 4"^)(1„3„2,-) . (4.58) 
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Following [55], we organise the infrared pole structure of the NNLO contributions renor- 
malised in the MS scheme in terms of the tree and renormalised one-loop amplitudes such 
that, 



Voles (4'^°) (1„ 3„ 2,-) + 4^^^) (1„ 3g, 2,-)) 



= 2 



r- ie) 



_,^r(l-2e) ((5, 



r(i-6) V 6 



K]l^^h2e)+H 



ma 



(2) 

mg 



+2/W(e)4^^°Hl„3„2,). 



(4.59) 



Here, 



(e) = AT (/«(e, .13) + /«(e, ^23)) - (e, S12 



+^^(4i!F(^'^13)+/S!F(e'^23)), 



with the individual if^^ defined in [32] and 



H 



(2) ^ 

4er(l-e) 
+ 



589 llvr^ 
432 ~ 72 



.rZ , 1. 41 TT^ 



3 tt" 



1 

7V2 



+ - 



19 TT 

18 ^36 



(4.60) 



54 24 / AT 27 



(4.61) 



We denote the finite contributions as, 



.;^mzte(4'"°) (1„ 3„ 2,-)) = Ar2^('x0),/mite ^ ^(2x0),/mite ^ 1 ^(2x0^,/imte 



iV2- 

4 (2x0), /mite .(2x0),/mjfe 



-L/\r2 aC^xO), finite j.j 



4 

N 



^ _ AT 1 (4.62) 



2 ^(lxl),/imte 



3,1 

+ArAr^4;^;^;^'"**^ + i^^ 

j_Ar2 /I (1x1), /inite 



1 ^(lxl),/imte 
7V2 

{lxl),finite , A^F .(lxl),/mtte 
jY ^3,AfF/iV 



(4.63) 



Explicit formulae for the finite remainders have been given in [18]. These are expressed 
in terms of one-dimensional and two-dimensional harmonic polylogarithms (HPLs and 
2dHPLs) [59,62], which are generalisations of the well-known Nielsen polylogarithms [63]. 
A numerical implementation, which is required for all practical applications, is available 
for HPLs and 2dHPLs [64]. 
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Finally, a finite NNLO contribution arises from the squared one-loop amplitude (4.29) 
for 7* ggg: 

A ......... ^^^^^^ 



Ct'\h, 2„ 3,) = N,,,(^^-n] ciX!f '''''' ■ 



5. Construction of the NLO subtraction term 
Three-jet production at the leading order is given by: 



(5.1) 



This leading order cross section defines the normalisation for all higher order corrections 
discussed in the following. 

At NLO, the tree-level four-parton processes 7* qqgg, 7* QQq'q' (non-identical 
quarks) and 7* — qqqq (identical quarks) yield three-jet final states. Only the two former 
processes require subtraction, since the third process is infrared finite. 

The four-parton real radiation contribution to the NLO cross section is 

dcr^LO = ^4 d$4(pi, ■■■,P4;q) 

I (i,i)GP(3,4) 

The antenna subtraction term is then constructed as: 



(5.2) 



d(TArLO = -^4 d$4(pi, ■■■,P4;q) 



(ij-)eP(3,4) 



4i2g,ig,jg) A^l g , (ji)„, (2i)-) J^ipi , , 



^3(1.,^., 2,-) Al{{li)^,jg, (2i)-) 4'\ni,Pj,P^i) 



2N 



El{lg,3,,Ag')Al{{13)^,{43)g,2g)4'^ip{s,^s,P2 



+El{2g, 3g,Ag') Alilg, (34)^, (24)-) (pi,m,m 

Integration of this subtraction term over the antenna phase spaces yields: 
dt^A^LO = ^3 (^) d$3(pi, • • • ,P3;q) 



(5.3) 



xAl{lg,3g,2g)4''\pi,P2,P3 



(5.4) 
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Together with the virtual one-loop contribution to 7* — qqg, 
\ [A\{lq, 3g, 2,-) + ^^(si23)A°(l„ 3„ 2,-)] 



^ iJa,, 3^, 25) +^^(si23)^0(lg, 3^,25-) +iV^ii(l„33,2g-)) , (5.5) 



N 



one obtains 



PoZes (da^^o) + Voles {AaJjLo) = > (5-6) 
thus yielding an infrared-finite result. 

6. Construction of the N"^ colour factor 

The N'^ colour factor receives contributions from five-parton tree-level 7* qqggg, four- 
parton one-loop 7* qqgg and tree-level two-loop 7* qqg. The multiple gluon emissions 
are colour-ordered, and the squared matrix elements do not contain interference amplitudes 
between different orderings. In the loop contributions to this colour factor, non-planar 
momentum arrangements are absent. 

6.1 Five-parton contribution 

At leading colour, the five parton contribution to three-jet final states arises from the colour- 
ordered emission of three gluons in 7* — > qqggg- The matrix element for one ordering is: 

= N,N^AUlg,ig,jg,kg,2g) . (6.l) 

The real radiation contribution to the cross section is obtained by averaging over all possible 
six orderings: 

df^^iVLO,JV2 = -^5iV^d$5(pi,...,p5;g)^ ^ A^{lq,ig,jg,kg,2g)J^^\pi,...,p5) 

■ (ij,fc)6P(3,4,5) 

[^5(1?' ^ff'is' ^9' %) + ^5(1?' ^ff'iff' ^ff' %)] "^3 ''(Pl) ■ ■ ■ jPs) ) 

{i,j,k)ePc{3A,5) 

(6.2) 

where the second expression is obtained by restricting the summation to the three cyclic 
permutations of the gluon momenta, while making the corresponding three non-cyclic per- 
mutations explicit. The cyclic form (6.2) is more appropriate for the construction of the 
real radiation subtraction term, since this form matches onto the full quark-gluon antenna 
functions of [32], which have a cyclic ambiguity in their momentum arrangements. 
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The real radiation subtraction term for this colour factor reads 




-ftihdg, kg) dl{\, (ij),, j ^3((lU%, S^),, 2,-) 4'\mj,P^uP2) 

-4{h, kg,jg) 4{mg. Wg, ig)^ Al{{lkj)g, (ijk)^, 2g) jf (pl^' , , P2 ) 

-(f^{lq,ig,jg)4{{^i)q,kg,{ji)g)^ A'i{{lik)g,{jki)g,2q) {p^k,Pjki,P2) 

+ (^DlaiU^ig,jg,kg) 

-dUh, kg,jg) 4{mg. ig. {jk)g)^ ^^(0^),, (7^),, 2,-) 4'\p^i,pJfk,P2) 

+ (^D%{2q, ig,jg, kg) " d°(2„ igjg) ( (2i) , , (ji),, ) 

(^., kg) 4{2q, (ij)g, (kj)g)^ Al{(2ij)g, i^g, l^") , , l ) 

+ i^l^{2q,ig,jg, kg) " /g^ (^„ , ) 4{2q, (kj)g, Jij)g) 

-4{'^g,kgJg)4{(2k)gM)g,ig)]A'3i(^)g,mg^^^ 
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+ \^Dl^{2g,igJg,kg) 

-4{2g,ig,jg) ( (2i) ^ , kg, (jij^j^ ( (2^^) ^ , (jki) ^ , Ig) {p^k , Pjj^i, Pi) 

-4(2^, kg,jg) 4{(2k)^, ig, (jk) ^^((2^1),, {jik) g , l,") {p^^ , pjfk , Pi) 

- (^Alil,, ig, kg, 2g) - ^^(l,, Z„ 2,- ) ( (H) ^ , kg, (2ijg) 

-A'i{lg,kg,2g)Al{(lk)^,ig, {2k)g)^ ( (Hfc) ^ , j„ (2^^ -) , P„ ) 

-^^3(l9,^<?>i5)rf3(2$> kg, {ji)g) AliiU)^, mk)g, (2k)g) jf\m,m)k,P^k) 

-^^0(2,-, kg,jg)4{l„ig, (Jk) g) A0((H),, ( W^),, (2^),-) jf\m,P^)M) 

fc„ j3)fi^(2,-, i„ (Tfc),) (c^)^, (2ij^-) 4^\prk,P^i,P^) 

-^d0(2,-, fc„ (ji)^) Aliilk)^, ({ji)k)g, (ij^-) 4\p{k,P^k,pii) 

+ ]^4{lq,ig,jg)4{jU)^,kg,(ji)g)Al{({^ 
+ l4i^,kg,jg)4m)^,ig,{jk)g)A'iiilj^ 

+l4i'^g,igJ9)dli(2i)g,kg, {ji)^) ("(jjfc)^, ({2{)k)^) {pi,p^k,P^k) 

+ ^4{2g,kg,jg)4{(2k)^,ig, (jl) , ) ( 1 ^ , ( , , (^2fe)i)-) jf (^1 , 

-lAO(l,,i,,2,-)dO((Hj^,A;,,j,)AO((a5^)^,(^^^^ 

+^d°(l,,fc3,i,)A|]((Tfe)^,i,,2,-)A0((a5i),,(7^^^ 

/c,, 2,-)dO((l^)^, ^^((05^),, (iij,, {2k)^) jf (^^pT^.p^fe) 
A;^, 2,-) (j^^ 
-^^0(l,,i,,2,-)d0((2i)-,/c,,j,)A0((lij^,(^^^ 

+ ^^3(2,-, kg,jg)Al{lg,ig, {2k)-) A^jU) ^, (jk) ^ , (^) {p{i , Pfk , P^i) 

-^^0(1,, /c,, 2,-)a!0((2fe)-, A0((T^^^ 
+^ciO(2,-,i5,j^)AO(l,,A;,,(20,0^3((l^),>(^^ 
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(6.3) 



6.2 Four-parton contribution 

The leading colour four-parton contribution comes from the one-loop correction to 7* 
qqgg, where the gluonic emissions are colour-ordered. It reads 

VA 



da 



NNLO,N'^ 



(6.4) 



(jj)e(3,4) 

The one-loop single unresolved subtraction term for this colour factor is 



da 



VS,1 

NNLO,N^ 



O d$4(;>i,...,P4;g) 



(jj)e{3,4) '- 
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Al{\,ig,3g,2q) jf'{pi, ■ ■ ■ ,Pi) 



1 



+ S 2 ^ 4(.U,h^J9) [^3((l^),,(i0,,2,-)+^^(sl234)^^((l^),,(J^),,2,-] 

^ (i,j)6(3,4) L 



+ 



\ {Vl{suj)+Vl{s^^^)) dl{\,tg,jg)Al{{li)^,{ji)^,2q) 



+ + \j^^{Sij) + \VI{S^^) - VliSUj)^ 4{l,,ig,jg)Al{{U)^, {ji)g,2g) 

+ 60 log ■^d'iil^igjg) A'iiiU)^, 2,-) 

Slij 



d°(l„i„i,)A0((H)^,(ji)^,2,-) 



jf\pu,Pji,P2) + (1^2) 



-- Yl Al{l,,ig.'^,)Al{{lt)^,j,,{2i)^) 
{2j)e(3,4) L 

+ (^^(512) -^§(.12^)) A0(l„z„2,)^g((Hj,, j„(2^,) 



Al{l,,^g,2q)Al{{l^)^,jg,{2^)^ 



Jz\pii,P2i,Pj) 



(6.5) 
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6.3 Three-parton contribution 

The thrcc-parton contribution consists of the two-loop three-parton matrix element to- 
gether with the integrated forms of the five-parton and four-parton subtraction terms, 



VS,1 

NNLO,N2 



+ ^^^3(^13) + ^^^3(^23) -^3(^12 



+ 77 (^?3°(^13) +2^^(523)) ^J(l„3„2,-) 



2e 



Vl{sx^) ((sia)-^ - (S123)-') +25^(S23) ((S23)-' - (S123)-') 



^0(l„3g,2,-)|da3, 
(6.6) 

(6.7) 



where we defined the three-parton normahsation factor 

— j d$3(pi,^?2,P3;g)4 (Pl'^'2,P3)- 

Combining the infrared poles of this expression with the two loop matrix element, we 
obtain the cancellation of all infrared poles in this colour factor, 

Voles {da%^j^Q^^^ + Voles (daJJ^^ j^^) + Voles (daJJ^^Q^^^) = . (6.8) 
7. Construction of the 7V° colour factor 

The contribution for the colour factor to three-jet final states is more involved than all 
other colour factors. It receives contributions from all partonic subproccsses: 7* — > qqggg 
and 7* QQQqg at tree-level, 7* qqgg and 7* — > qqqq at one loop and 7* qqg at two 
loops. All contributions contain a mixture of colour ordered and non-ordered emissions. 

7.1 Five-parton contribution 

The subleading colour N'^ contribution of five-parton final states to three jet final states is 
do"^iVLO,Aro = ^5 d$5(Pi, ■ ■ ■ ,P5; 9) 

iO/l Q /I K O \ \ ^ a0( 



1 

3! 



^5(l?5 3fl>4g,,5g,2g) ^ A'^{lg,ig,jg,kg,2q) 

(ij,fc)GP(3,4,5) 



(-'-95 '^95 ^95 ^q) "I" ^5 (^95 '^9' ^qi ^9' Iq) ^5 (-'-q' '^9' ^gi ^q) 
= Ar5d$5(pi,...,p5;g) 



■4^^(Pi, • • • ,^5) 



(i,i)GP(3,4) 
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'^^bi^qj ^qj ^gj + C^{2q,Aq, 3q, 5g, Ig) — C^{lq,3q, Ag, 5g,2q) 



(7.1) 



where the symmetry factor in front of is due to the inherent indistinguishabihty of 
gluons. In A^, gluon (5g) is effectively photon-hke. It does not participate in any three- 
gluon or four-gluon vertices, and there are no simple collinear limits as («)g||(5)g and 
(J%ll(5)p. 

The real radiation subtraction term for this colour factor is: 



^'^NNLO,N0 



N5N^d^5{pi,---,P5;q) 



'I E A^s{lg,5g,2g)Al{{l5)g,ig,jg, {25)^)4''^ {pi5,Pi,Pj,P25) 

Al{{li) (ji) 5g,2q)4^\p{i,p]'i,P5,P2) 



+ 4{2, j, i)A'i{lg,5g,iij)g, (2 j) -) ^ {Pl , P5 , Pij , P2j ) 

1 



+ ^ Yl A^3{U^ig,2g)Al{{li)^,jg,kg,{2i)^)4'^\pii,Pj,Pk,P2i) 

■ (i,i,fc)GPc(3,4,5) 

-A'i{lg,5g,3g) [c^Olb)^, (35) ^ , ^ , 2,-) + C4°(2,-,4,-, (35)^, (Ts)^)] 4^\p^5,P^5,P4,P2) 
-Al{2g,5g,4g) [c'4^1,,3,,(45)^-,('25)^-) + C4°(('25)^-,(45)^-,3,,^ ^^(^1,^3,^^5,^15) 

4 E ( ^4(1?, igJg, 2,-) - ^^^(1,, igjg) ^^((H),, {ji)g, 2,") 



(i,i)6(3,4) 



-d°(2,-,j„z,)^0(l„(zj%,(2j),-)j Al{{lij)^,5g,{2ji)g)4'\p^j,p,,p^) 
-\ E ( ^4(1., ^s, 5„ 2,-) - ^°(1„ z„ 2,-) ^°((n),, 5„ (20,-) 



(j,j)e(3,4) 



-A°(l„5„2,-)A0((15),,i„(25),) 1 Ali{1^5)^,jg,{2^5)^4'\m5,PJ,mf5) 



+ - i°(l„ 3„ 4g, 2g) - A'iilg, 3g, 2 g) ^^((l3)^, 4g , (23), 



-A0(1,,43,2,-)^0((14),,3„(24),-) 1 A^((134)„, 5^, (234)„-) jf (pir4,^^5,P^) 



[C4°(l„ 3„ 4,-, 2,-) + C4°(2,-, 4g, 3g, 1,)] Al{{13A) 5g, (234) ) jf {p^a,P^m) 



+^ E I ^3(l9,%>ip)^3((l0,,5,,25-)yl0(((H)5),,(jz)^,(25),-) jf (^5,P7i,P^^ 



(i,i)G(3,4) 



+ 4{2g,jg,lg)Al{lg,bg, ( 2^),-) ^3 ( ( 1 5 ) , , (U^, ((2j)5),-) (^^5 , P^' , ^^^5 



(3), 
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7.2 Four-parton contribution 

The four-parton contribution comes from the subleading colour one-loop correction to 
7* ~^ QQ99j where the gluonic emissions are colour-ordered, from the leading colour one- 
loop correction to 7* <l<19g, where the gluonic emissions are not colour-ordered and from 
the one-loop correction to the identical-flavour process 7* qqqq- It reads: 



da 



v,i 



iV7VLO,iV0 =^4iV° (^) d$4(pi,...,P4;g) 



2^ ^ i^q^h^jg^'^q) ^4' (Igi ^gi ^sj ^q) -|- (1^, 3^, 4^, 2^) j 

^ (i,i)e(3,4) ^ 

-Cl''^(l„3„4,-,2,-) - Cl''^(2,-,4,-,3„ Ji'^(pl,P2,P3,^>4), (7.3) 
The one-loop single unresolved subtraction term for this colour factor is 



(ij)e(3,4) 
+Al{si2)Al{lq,ig,jg,2q) 



\2tt 



4^^(Pl,P2,^'3,P4) 



-^^0(si2)i^(lg,3g,43,2,-)4')(pi,P2,P3,P4) 

+ {Alis^^) + AI{s2a)) (C0(l„3„45,25) + C0(2,-,4,-,3„l,)) 4*\pi,P2,P^,pd 
1 
2 



dl{}q,ig.3g) i^((H),,(ii),,2,-)+^i(si234)A°((lz)^,(iz)^,2,-) 



(«J)e(3,4) 
■^3 {Pu^PjhP2) + (1^2) 



- Alilg, ig, 2g) Al{{li)^,jg, (2z)„-) jf\piuPj,P2i) 



-Al{lq,ig,2^) A\{{li)^,jg,{%)-) + Al{si2M)Al{{lt)^,jg,{2j)-) 4\m,Pj,m 



-A\{lg, ig, 2g) ( ( li) , , j„ {2i) -) (pTi , P,' , P5 



d°(l,,z,,j,).4°(s^J^°((li)^/(ji),,25)jf (pr^,g'i,P2) + (1 - 2) 



(li)2 



[^|](si2) - ^^(^120] Alii,, z„ 2,-) Al{{lt)^,jg, {2^)^ J^i\m,P,.m) 



Vl{s~) + Vl{s~) ^0(l„z„2,-)A0((lz),,jg,(2z)^-)jf (^i,i?,-,^i) 



Al{lg,ig,2q)Ai{{li)^,3g,{2^)^Jf\n^,p,,n. 



+ 
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i,, 2g) Al{{u)^,jg, {2i)^ 4\p{i,Pj,m 



-bolog^Al{l„i„2^)Al{ili) jg,{2i) )jPi^,,p„^^\ (7.4) 

Sl2i 



7.3 Three-parton contribution 

The three parton contribution to the colour factor receives contributions from the three- 
parton virtual two-loop correction and the integrated five-parton tree-level and four-parton 
one-loop subtraction terms, which read 

^^nnlo,no ^"nnlo,no ~ -^^ 



-AI{S,2) AI{S,2) + AI{S12) + A\{S,2) 



-\ (^3(^13) + 'Dl{s2^)) A\{1^, 3„ 2,-) - Alisi^) A\{1„ 3„ 2,-) 

-^^^(^12) ((512)-^ - (5123)-^) ^§(l„3„2,-)|da3 . (7.5) 

Combining the infrared poles of this expression with the two loop matrix element, we 
obtain the cancellation of all infrared poles in this colour factor, 

Voles (dcj;^^^o^^o) + Voles (daJJJJ^Q^^o) + Voles (daJ^^'J^^^^^o) = . (7.6) 
8. Construction of the colour factor 

The I/N"^ colour factor receives contributions from five-parton tree-level 7* qqggg and 
7* ~^ QQQQ9j four-parton one-loop 7* qqgg and 7* qqqq as well as tree-level two-loop 
7* qqg. The gluon emissions are all photon-like, not containing any gluon self-coupling. 
The four-quark processes contribute through the identical-quark-only terms. 

The construction of the subtraction terms for this colour factor was discussed in detail 
in [32,65]. 

8.1 Five-parton contribution 

Two different five-parton final states contribute at 1/7V^ to three-jet final states at NNLO: 
7* — qqggg and 7* — > qqqqg with identical quarks. 

The NNLO radiation term appropriate for the three jet final state is given by 



^'^NNLO,l/N^ = ^ d^>5(pi, . . . ,P5; g) 



^ ^ti^q^ ^3' ^g^^gi 2«) + 2(55(1^, 2q, 3q,, 4^, 5g) 
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where the symmetry factor in front of A'^ is due to the inherent indistinguishabihty of 
gluons. The factor 2 in front of arises from the fact that two different momentum 
arrangements contribute to the squared matrix element (4.21). If the quarks and antiquarks 
are not distinguished by the jet functions, these contribute equally. 
The real radiation subtraction term for this colour factor is: 

<i<^NNLO,l/N^ = ^d$5(pi, . . . ,P5;q) 

1^ H ^3(l9>^9,2,-)i^((H)^,i5,A;5,(2ij-) jf)(pri,Pj,Pfe^ 



+ {A'i{lg,ig,jg,2g) - ^^(l„i„2,-) ^^((iz),, J„ {2i) ^) 

-A^s{lg,jg,2g) A^,{(lj) ^, ig , (2j) -) ) ^^((lii)^, A;,, (2ii) -) jf (pl^, Pfe, pj^) 



+2 



^0(1,, 5g,2g) C0((15)g, 3„ 4g, {25) 4''\m5,P3,P4,m5) 

+A°{lg, 5g, %) C2{{l5)g, 3g, (45)^-, 2g) jt\p^5,P2,P3,P^5) 
+Al{3„ 5g, 2g) C2{lg, (35)^, 4g, (25)^) jt\m,Pi,PA,m5) 
+Al{3„ 5g, 4g) C0(1„ (35)^^, (45)^-, 2,-) (^^5,^1,^2,^) 

-Alilg, 5g, 3g) C ^ ( (iS) ^ , (35)^, 4g , 2g) (^„ ^2 , ^4 , ) 

-^^(2^, 5^, 4g) C2{lg, 3g, (45)^-, (25)^) 4''\m,PuP3,m) 
+C2(.lg, 3g, 4g, 2g) AI{(134) ^ , 5g, (234)-) Jf (pIl4,P5,P^4) 



(8.2) 



The sum in the first contribution runs only over the three cyclic permutations of the gluon 
momenta to prevent double counting of identical configurations obtained by interchange of 
j and k. 

8.2 Four-parton contribution 

At one-loop, there are two contributions to the colour suppressed contribution proportional 
to one from the four quark final state and one from the two quark-two gluon final 

state: 

d^SiO.i/JV^ = ^ (^) d$4(pi, . . . ,P4; g) 

I ^ Al'\lg, 3g, 4g, 2g) + 2 C^ilg, 3g, 4g, 2g)^ J^^^ (^1, . . . , ^4 ) , (8.3) 

where the origin of the symmetry factors is as in the real radiation five-parton contributions 
of the previous section. 

The corresponding subtraction term is: 

^''InIomn^ = ^ (^) ^Mpu ■■■,P4;q) 
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I ■ i,ieP(3,4) L 

+ (Al{lg,ig,2g) ( (H) ^ , , ^i) -) + 1234 ) ( , , J, ^ 

-2 [At{s^2) + AUsu) + At{s2s) + A^{ss4) - Atisis) - At{s24 

XC^{lg,3g,4g,2g)4''\pi,PS,P4,P2) } ■ 



(8.4) 



8.3 Three-parton contribution 

The three-parton contribution consists of the two-loop three-parton matrix element to- 
gether with the integrated forms of the five-parton and four-parton subtraction terms, 



da 



VS,1 



1 



NNLO,l/N^ ^"NNL0,1/N'^ ~ ^2 



2^4(^12) + 2C4°(512) + AI{S12) Al{lq, 3^, 2,-) 



+^°(si3)ikl„3„2,-)lda3 



(8.5) 



Combining the infrared poles of this expression with the two loop matrix element, we 
obtain the cancellation of all infrared poles in this colour factor, 

Voles (d(7;^jvLO,i/iv2) + 'Poles {da]fj^\Q y^^ + Voles (da]f^^Q^^^^^ = . (8.6) 
9. Construction of the Np N colour factor 

The colour factor Np N receives contribution from the colour-ordered five-parton tree-level 
process 7* qqq'q'g-, the four-parton one-loop processes 7* qqq'q' and 7* qqgg at 
leading colour and from the two-loop three-parton process 7* qqg. 

9.1 Five-parton contribution 

The NNLO radiation term appropriate for the three jet final state is given by 
^^NNLO,NfN = 

N^NFNd^^{pi,...,p^;q) 
5°'»(1„ 5g, 3,., 2,-) + (1„ 3,., 5^, 2,-)] (pi, . . . ,p^)i^.l) 

The two terms represent the two colour orderings of the leading colour amplitude for 
this process. Since the leading coloTir gg'q'g'-antenna subtraction terms allows q to represent 
either a quark or an antiquark, both colour orderings are mixed together. Therefore, it is 
not possible to construct a subtraction term for an individual contribution, but only for 
their sum. 
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The subtraction term for this contribution is 

'i<^ffNLO,NpN = ^5 NpNd^sipi, ...,P5;q) 

X {AO(l,,55,450SO((T5)g,3,^(4r)^-,,25) jf)(^5,P2,P3,1^5) 
+Al{3g,,5g, 2g) b2{1„ (35)g„4,-,, (25)^-) 4''\pi,p^5,P^5,P4) 

+Gl{5g, 3,,, 4,-0^^(1,, (U)g, 2g) jt\pi,P2,m4,P^A) 

-Gl{5g, 3,,, 4,-0 4{1„ (34)^, (54)^)) Al{(m\, (543)^, 2,-) jf (^^14,^2,^13) 
+ (El,ilg, 3g,,4g',5g) - Al{l„ 5„ 4,-,) El{il5)^, 3,,, (45)^-,) 

-G°(5,,3,^4,-0d^(l5,(54)^,(34)^)) ^^(0^),, (S^^, 2,-) (pir4,^'2,pS5) 

-G°(5„ 3,,, 4,-0 d0(2,-, (43)^, (53)^)) (534)^, (243)^04'^(Pi,;>^3,;>^4) 

+ (^^0,(2,-, 4,-,, 3,,, 5,) - Ali3g,,5g, 2g) El{(2^)^Aq', (35) ,0 

-Gl{bg, 3,,, 4,-0 dg(2,-, (53),, (43),)) ^^(l,, (435),, (253),0 jf (pi,Kr3,pir5)} ■ 

(9.2) 

9.2 Four-parton contribution 

The four parton contribution to the NpN colour factor reads: 

'^'^NNLO,NpN = ^4 NpN d$4(pi, ■■■,PA]q) 

' \ XI (^4'''(l?'V'J9''2q) -Fvl4'''(lg,ig,ig,2^)H 4'^^(pi,...,]34). 

. (i,j)e(3,4) J 

(9.3) 

The average over the permutations of the momenta (3) and (4) has to be made in both 
contributions to this colour factor. In the one-loop correction to the 7* — qq'q'q final 
state -B4'", the secondary quark-antiquark pair has to be symmetrised, since the quark- 
gluon antenna functions used in the one-loop subtraction terms do not distinguish quarks 
and antiquarks. The summation over the two colour orderings of the one-loop correction 
to the 7* qggq final state A^'^ must be kept since the one-loop subtraction functions 
appropriate to this term contain both orderings because of their cyclicity. 
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The corresponding subtraction term is: 



'^^NNLO,NirN = NaNfN j d$4(pi, • • • ,P4; q) 

I fj,7)e(3,4) 



(«,i)e(3,4) 

(ij)e(3,4) 



+ <! o (^3°(1.,3,.,%) fA^((13)^,(43)^,2,-)+^i(si234)^g((13),,(43)^,2,, 



3,,,V) ^1]((13)^, (43)^, 2q) 
+\ (253°(si34) + 253°(^2{45))) ^3°(1., 3,., 4,-,) ^°((13)^, (43)^, 2,-) 
+ (^^(sia) + ^^(si4) - P§(5i34)) £^3°(1„ 3,., 4,-0 A0((13)^, (43)^, 2,-) 

+6olog^^3^1g,3g^4g-0^^(CL3) (43) 2,-)") +(1^2) 
Sl34 y a y 

+ I i f^3(l.> 3„ 4,)i^((13)^, (43),, 2,-) + Di(l„ 3„ 4,) ^^((Ts),,, (43),, 2,-) 



+6o,Flog^D30(l„3„4,) A0((13),,(43),,2,-) ) jf (pr3,pr3,P2) + (1^2)^. 



+2^0(534) Dl{l,, 3„ 4,) A0((13),, (43),, 2,-) 



Sl34 



(9.4) 



9.3 Three-parton contribution 

The three parton contribution to the NpN colour factor contains the three-parton vir- 
tual two-loop correction and the integrated five-parton tree-level and four-parton one-loop 
subtraction terms, which read 



^^nnlo,NfN + 



vs,i 

NNLOMfN 

1 



+ 2^1(^23) 



^0(1„3„2,-) + - {£lisn) + £lis23)) ^^(l5,3„2,-) 



+ 2(^3(si3)+2?3°(523)) ii(l,,3„2^ 



bo,F 

' 2e 



^53°(^13) ((^13)-' - {S123D + ^3(^23) ((^23)-^ - M^') 



^3(lgi 3g, 2,) 
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2e 



SiM ((513)-^ - (5123)-^) +^3^(^23) ((^23)-^ - (^123)-^ 



A0(l„3^,2,-)|da3. 

(9.5) 



Combining the infrared poles of this expression with the two loop matrix element, we 
obtain the cancellation of all infrared poles in this colour factor, 

Voles (da^jvLO.AT^iv) + 'Poles (da^^^^ j^^j^^ + Voles (da^l^^^^^^^ = . (9.6) 
10. Construction of the Np/N colour factor 

The Np/N colour factor receives contributions from five-parton tree-level 7* qqq'q'g, 
four-parton one-loop 7* qqq'q' and 7* qqgg at subleading colour as well as three- 
parton two-loop 7* qqg. The gluon emissions are all photon-like. 

This colour factor is part of the QED-type corrections. We described the construction 
of the subtraction terms for this colour factor previously in [66]. 

10.1 Five-parton contribution 

The NNLO radiation term appropriate for the three jet final state is given by 

d'^NNLO,Np/N = ^'5^d$5(pl,...,p5;5) S°''^(lg, 5g, 2,; 3g/, 4gO 

+B^/{lg, 2g; 3g,,5g, 4,-0 - 2 B°'^(l„ 2,-; 3,,, 4,-,; 5,)] 4'\pi, . . . ,^5) 

(ij)e(3,4) 

+B^/{lg,2g;ig,,5g,jgr) " 2 s"'" ( 1„ 2,-; z,, , j,-, ; 5g )] ^ (^1 , . . . ,^5) , 

(10.1) 

where the symmetrization over the momenta of the secondary quark-antiquark pair ex- 
ploits the fact that the jet algorithm does not distinguish quarks and antiquarks. This 

symmetrisation reduces the number of non-vanishing unresolved limits considerably, since 
the interference term in S^'*^ is odd under this interchange. As a resTilt, the unresolved 
structure of the symmetrised i^^'^ equals the unresolved structure of B^'^ + B^''^. 



S 



The subtraction term reads: 

(pi,...,P5;g)^ J2 \ 

-A^lq, Bg, 2g) b2{{i5)^, (25)^, jq') jtHp^5,P^5,Pi,Pj) 
-Al{ig>,5g,jq>) B2{lg,2q,{i5)g,,{j5)g,) J^^^ (pi , , ^ , P^s ) 

{El{lg,iq,Jg,)Ali(li)g,{ij)g,B9^'^g) 4'\p^J^P2,Prj,P5) + (1 ^ 2)} 

-( 50(1,, v,J,-',2,-) - i {El{lg,ig,Jq,)Al{{li)^,iji)g,'^g) + (1 2)} 
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xA0((li5)_, (i5i)„, 2g) 4\^5,V2,vJii) + (1^2) 



-ii £;3°(i,, V, i^O ^slOiT., 2,-) A0((ai)5 



+(1 ^ 2) 



(10.2) 



10.2 Four-parton contribution 

The four parton contribution to the Np/N colour factor reads: 



{-I E (sl''(l„V>J.-'>2,-) + 2Cl'^(l„z„i,-,2g 

+^4 ""(Ig, ^9,^3, 2^)) I ji^Hpi,...,P4)- 



(10.3) 



Like in the Np N colour factor, the expression is symmetrised over the momenta (3) and 
(4) to remove terms which are antisymmetric under charge conjugation, and can not be 
accounted for properly by the quark-gluon antenna functions. 
The corresponding subtraction term is: 

Np /as\ 



+ ^ [£lisi3) + Slisi^) + 6lis23) + SI{S24)] i°(l„ 3„ 4„ 2g) jf (pi, . . . ,P4) 
EliU, 3,., 4,-0 \Ali{13)^, (43)^, 2g) + ^i(.i234)^^(('T3),, (43),, 2g) 



+-^3(^(13)2) ^3(1?, 3,., 4,-0 A0((13),, (43),, 2,-) 



Elilg, 3,,, 4,-0 + ^°(s34)i?3'(l9> 3g', 4,-0 ^E((13),, (43)„, 2,0 ) (m3,pr3,^'2) 



+(1 ^ 2) 



- ^ (A0(l„z„2,-)ii((li)^,j„(2i)-) + 

(ij)e(3,4) ^ 
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+h,F\og^Al{l,,ig,2^)Al{^)^,jg,'{^^ (10.4) 



10.3 Three-parton contribution 

The three parton contribution to the Np/N colour factor consists of the three-parton 
virtual two-loop correction and the integrated five-parton tree-level and four-parton one- 
loop subtraction terms, which read 

A S I J VS,l _ 

'^^NNLO,Np/N ^^NNLO,Nf/N " 

+A{si2) + ^4'(^i3) + ^4'(^23)1 Al{l^,2>g,2q) - \ (^3°(^i3) +^3(^23)) ii(l„3„2,-) 



2 



-^0(5i2)ii(l„3„2,-) -hLAl{s^2) {{S12)-' - (5123)-^ ^°(l„33,2,-)|da3 . (10.5) 

Taking the infrared pole part of this expression, we obtain cancellation of all infrared 
poles in this channel: 

Voles {daf^NLO,Np/N) + ^o/es (daJ^^LO.iV^/iv) + ^^^^^ {^''^nIlo,Np/n) = • (10-6) 
11. Construction of the Np colour factor 

The N"^ colour factor receives contributions only from the four-parton one-loop process 
7* ~^ QQQ'q' aiid from the three-parton two-loop process 7* qqg. 

This colour factor is also part of the QED-type corrections, described previously in [66] . 

11.1 Four-parton contribution 

The four-parton one-loop contribution to this colour factor is 

'^''NNLO,Nj. = ^^^F (^) d$4(pi,...,P4;g)i?l''(lg,3g',V,2,-) jf)(pi,...,p4).(ll.l) 

This contribution is free of explicit infrared poles (as can be inferred from the absence of 
a fivc-parton contribution to this colour structure). 

The subtraction term corresponding to this contribution is 



J vs,i 

^^NNLO,Nj. 



N^Nl (^) d$4(Pi,...,P4;g)^| 



El{l,, 3,,, %) + 6o,F log ^ El(l,, 3,,, %0 ^°((13)^, (43)^, 2,-) 

S134 
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•4^^(P13,P43,P2) 



+ 



El{2g,3,,Aq') + bo,F log^El{2g,3,,,4^,)) (43)„, (23) 



S234 



(11.2) 



Although and ^3 contain cxpUcit infrared poles, these cancel in their sum, as can be 
seen from (5.16) and (6.32) of [32]. da"]^'^^^ ^2 is therefore free of explicit infrared poles. 

11.2 Three-parton contribution 

The three parton contribution to the Np colour factor consists of the three-parton virtual 
two-loop correction and the integrated four-parton one-loop subtraction term, which reads 



da 



VS,1 

NNLO,Nl 



4'(si3)+4'(^23)) Al{l„3g,2g) + {Slisn) + Slis23)) ij(l„3^,2,-) 
+^ [£l{s,s) {{s^s)-^ - (.123)-^) +£I{S23) {{S23)-' - {S123)-')] Al{l„3„2,) 



dcra . 
(11.3) 



Combining the infrared poles of this expression with the two loop matrix element, we 
obtain the cancellation of all infrared poles in this colour factor, 



Voles ^dcj 



IfNloXi) (d^^LCivi 1=0- 



(11.4) 



12. Construction of the Nf,^ colour factor 

The Np^-y colour factor comes from the interference of amplitudes in which the external 
vector boson couples to different quark lines. It receives contributions from five-parton tree- 
level 7* — >■ qqq'q'g, four-parton one-loop 7* qqq'q' and 7* — > qqgg as well as three-parton 
two-loop 7* qqg and 7* ggg. This colour factor is absent in three-jet production at 
NLO and four-jet production at LO because of Furry's theorem [5], and its contribution 
to four-jet production at NLO is numerically tiny [52]. The numerical magnitude of this 
term in the two-loop corrections to the three-parton channel is equally very small [18,53]. 
A detailed discussion of this colour factor, and of the effects leading to its numerical 
suppression is contained in [52]. All partonic channels contributing to this colour factor 
are individually finite. 

In this section, we document this colour factor for completeness. Given that is numer- 
ical impact can be safely expected to be negligible, we refrain from its implementation. 
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12.1 Five-parton contribution 

The NNLO radiation term appropriate for the three jet final state is given by 

d<^NNLO,NF,-, = ^5 Nf,^ d$5(pi, ■■■,P5;q) 

x[-N 5g,4gr,3g,,2g) + B^^'\l g , Aq' \' , ^ , 2g-) - 45,; 3g,, 25, 5^ 

+^ (4°'=(l5,5,,2,-;3,^%0 +4°'"(l,,2,-;3,^5,,4,-0 +4°''=(l5,2,-;3,,,%,;5,)) 
xjf (pi,...,p5) • (12.1) 

Once symmetrised over the quark and antiquark momenta, this term can be integrated 
safely without the need for an infrared subtraction. It is free from infrared singularities 
associated with gluon 5g unresolved, since the corresponding four-parton tree-level term 
vanishes after symmetrisation over the quark and antiquark momenta. Double unresolved 
singularities can not appear since there is no tree-level three-parton process proportional 
to Np^j- 

12.2 Four-parton contribution 

The four parton contribution to the Np^j colour factor reads: 



+ NAl'\lg,ig,jg,2g) - ^ A^/ilg , 3 g , , ^g) ^ (R , . . . , ^4 ) • 

(12.2) 

Terms which vanish under symmetrisation of the quark and antiquark momenta, arising 
from in (4.46) have been omitted here. After this symmetrisation, all explicit infrared 
poles present in individual terms in the above expression cancel. Moreover, (12.2) is finite 
in all single unresolved limits, such that no antenna subtraction is needed. 

12.3 Three-parton contribution 

The three-parton contribution to the Np^j colour factor consists of the three-parton virtual 
two-loop correction to 7* qqg [18] and the one-loop squared correction to 7* ggg [53] . 
Both are individually finite, and were shown to be numerically tiny. Since no subtrac- 
tions were carried out in the five-parton and four-parton channels, there are no integrated 
subtraction terms in the three-parton channel. 

13. Numerical implementation 

Using the matrix elements and antenna subtraction terms derived in the previous sec- 
tions, NNLO corrections to any infrared-safe three-jet observable in e^e~ annihilation (jet 
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Figure 7 : Structure of the EERAD3 parton-level Monte Carlo event generator programme. 



cross section, event shape variable) can be computed numerically. Wc implemented this 
numerical evaluation into a parton-level event generator program, which we name EERAD3. 

This program is based on the program EERAD2 [39], which computes four-jet production 
at NLO. EERAD2 contained already the five-parton and four-parton matrix elements relevant 
here, as well as the NLO- type subtraction terms dcr^^^Q and dcr]^^^^. 

The implementation contains three channels, classified by their partonic multiplicity: 



in the five-parton channel, we integrate 



NNLO 



• in 



the four-parton channel, we integrate 



da 



v,i 

NNLO 



da 



VS,l 
NNLO 



in the three-parton channel, we integrate 



da 



iVATLO + '^^NNLO + d<7 



VS,1 
NNLO 



(13.1) 



(13.2) 



(13.3) 



The numerical integration over these channels is carried out by Monte Carlo methods using 
the VEGAS [67] implementation. The structure of the programme is displayed in Figure 7. 

The phase space in the four-parton and five-parton channel is decomposed into wedges 
which are constructed such that two of the invariants are smaller than any of the other 
invariants. This decomposition allows an optimal generation of phase space points in 
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the unresolved limits. The full four-parton phase space is obtained by summing (a) 12 
wedges with {sij,Sik) smallest, plus (b) 3 wedges with {sij,Ski) smallest. To obtain the 
full five-parton phase space, we sum (a) 30 wedges with (sij,Sik) smallest, and (b) 15 
wedges with {sij,Ski) smallest. The phase space integration in either channel is carried 
out by integrating only over a single wedge of type (a) and a single wedge of type (b), 
while summing the integrands appropriate to all wedges of the given type. In doing this 
summation, we combine (in the exact unresolved limits) phase space points which are 
related to each other by a rotation of the system of unresolved partons, thereby largely 
cancelling the angular-dependent terms. In all colour factors containing angular-dependent 
terms, the combination of phase space wedges yields a substantial improvement of the 
numerical stability of the results. 

It was already demonstrated above that the integrands in the four-parton and three- 
parton channel are free of explicit infrared poles. In the five-parton and four-parton chan- 
nel, we tested the proper implementation of the subtraction by generating trajectories 
of phase space points approaching a given single or double unresolved limit using the 
RAMBO [68] phase space generator. Along these trajectories, we observe that the antenna 
subtraction terms converge towards the physical matrix elements, and that the cancella- 
tions among individual contributions to the subtraction terms take place as expected in 
the antenna subtraction method. 

Moreover, we checked the correctness of the subtraction by introducing a lower cut 
(slicing parameter) yo on all phase space variables, and observing that our results are 
independent of this cut (provided it is chosen small enough). This behaviour indicates that 
the subtraction terms ensure that the contribution of potentially singular regions of the 
final state phase space does not contribute to the numerical integrals, but is accounted for 
analytically. 

14. Thrust distribution as an example 

To illustrate the implementation and to study the numerical impact of the individual NNLO 
contributions, we consider the thrust distribution. We already reported the NNLO results 
on this event shape distribution in a previous paper [33], where the phenomenological 
implications are discussed in detail. 

The thrust variable for a hadronic final state in e+e~ annihilation is defined as [69] 



wheie Pi denotes the three-momentum of particle i, with the sum running over all particles. 
The unit vector n is varied to find the thrust direction Ht which maximises the expression 
in parentheses on the right hand side. 

It can be seen that a two-particle final state has fixed T = 1, consequently the thrust 
distribution receives its first non-trivial contribution from three-particle final states, which, 
at order ct^, correspond to three-parton final states. Therefore, both theoretically and 
experimentally, the thrust distribution is closely related to three-jet production. 




(14.1) 
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A study of the phenomenological implications of the NNLO corrections to the thrust 
distribution was presented in [33], illustrating that the NNLO corrections amount to about 
15% of the total result over the experimentally relevant range 0.02 < 1 — T < 0.3, and that 
inclusion of these corrections results in a considerable stabilization of the renormalisation 
scale dependence of the theoretical prediction. In the present context, we use the thrust 
distribution only as an example to illustrate certain features of our calculation. 

The three-jet rate and event shapes related to it can be expressed in perturbative QCD 
by dimensionless coefficients. These coefficients depend, for non-singlet QCD corrections, 
only on the jet resolution parameter (respectively on the event shape variable). Typically, 
one denotes these coefficients by A, B,C,... at LO, NLO, NNLO, etc. 

The perturbative expansion of thrust distribution up to NNLO for renormalisation 
scale /J? = s and ag = as{s) is then given by 



dT ~ V2^J dT ^ V2^/ dr^V2^J dT ' 



(14.2) 



Here we define the effective coefficients in terms of the perturbatively calculated coefficients 
A, B and C, which are all normalised to the tree-level cross section 



Ana ,^ 9 



(14.3) 



for e^e — > qq. Using 



=<.„(! + \C, (g) + (g)^ + Oid) ) , (14.4) 

with {Cf = {N^ - l)/(2iV), Ca = N,Tr = 1/2 for iV = 3 colours and Np hght quark 
flavours) 



..4 



123 



-C'f + CfCa - 44C3 + CfTrNf (-22 + I6C3) 



(14.5) 



we obtain: 



A = A 
B = B 



C = C- -Cf B + -Ci -K2] A. 



9 



(14.6) 



These coefficients depend only on the jet resolution parameter or the event shape variable 
under consideration, and are independent of electroweak couplings, centre-of-mass energy 
and renormalisation scale. 

The above coefficients include only QCD corrections with non-singlet quark couplings. 

At 0{af), these amount to the full corrections, while the 0{ag) corrections also receive a 
singlet contribution. As discussed above, this singlet contribution arises from the interfer- 
ence of diagrams where the external gauge boson couples to different quark lines. In four-jet 
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Figure 8: CoefBcients of the leading order and next-to-leading order contributions to the thrust 
distributions. 
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Figure 9: Coefficient of the next-to-next-to-leading order contribution to the thrust distribution. 
Solid: corrected for large-angle soft terms; dotted: original result. 

observables at ©(a^), these singlet contributions were found to be extremely small [52]. 
Also, the singlet contribution from three-gluon final states to three-jet observables was 
found to be negligible [53]. 

We determine A,B,C from the perturbative contributions to the differential cross 
section, normalised to the tree-level hadronic cross section: 

dT ~ (70 dT ' dT ~ (70 dr ' dT ~ ao dT ' ^ ' ' 

For the determination of the non-singlet coefficients, it is sufficient to consider ao for 
pure photon exchange, since any electroweak coupling constant cancels out in the above 
ratio. The LO and NLO coefficients A{T) and B{T) were computed in the literature long 
ago [5-8,10]. They are displayed for comparison in Figure 8. 

The total NNLO coefficient C(T) is displayed in Figure 9. The six different colour 
factor contributions to it are shown in Figure 10. It can be seen that the numerically 
dominant contributions come from the N'^ and Np N colour factors. The contributions of 
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1-T 1-T 

Figure 10: Different colour factor contributions to NNLO coefficient of tlie thrust distribution. In 
N'^ and colour factors: solid: corrected for large-angle soft terms; dotted: original result. 



these two colour factors are of opposite sign, with N"^ being of larger absolute magnitude, 
thus resulting in a total positive result. Contributions at the 10% level of the total come 
from the N'j^ and colour factors, Np/N amounts to about 5%, while the most subleading 
I/N"^ colour factor is below 1%. 

To illustrate the independence of our results on yo, we display the different colour factor 
contributions to C(T) as function of — ln(l— T) in Figure 11 for different values of the phase 
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Figure 11: Dependence on phase space cut yo in different colour factors. We see that the results 
are independent of yo for — ln(l — T) < 4, but, as explained in the text, differ at larger values of 
-ln(l - T) 



space cut yo = 10~^, 10~^, 10~^. By rescaling all phase space invariants to the total centre- 
of-mass energy squared, yo becomes dimcnsionless. Since the value of (1— T) determines the 
typical scale of the smallest resolved invariant, one must require yo to be several orders of 
magnitude smaller than (1— T) for the cancellation between matrix element and subtraction 
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term to be accurate. Figure 11 shows very clearly that over the phenomenologically relevant 
range, i.e. 0.02 < 1 — T or equivalently, — ln(l — T) < 4, our results do not depend on j/q- 
As (1 — r) approaches yo (starting at about (l — T) ^ C'(lOOO)yo)) the calculation becomes 
unreliable as expected. This behaviour can be understood to arise from the fact that the 
subtraction terms converge to the full matrix element only once all unresolved invariants 
are much smaller than any of the resolved invariants. 

The numerical convergence of our calculation deteriorates for lower values of yo for two 
reasons. 

• the absolute magnitude of matrix elements and subtraction terms increases for de- 
creasing yo both in the five-parton and four-parton channel. Consequently, numerical 
cancellations between matrix elements and subtraction terms happen over larger or- 
ders of magnitude, thereby enhancing numerical rounding errors. 

• the four-parton one-loop matrix elements start themselves to become numerically 
unstable because of the presence of inverse Gram determinants, which can become 
singular inside the integration region. 

Therefore, for all phenomenological applications to the thrust distribution [33], we 
choose yo = 10^^. For applications to other event shapes, one expects a similar behaviour, 
and one must first determine the value of yo required for reliable predictions in the phe- 
nomenologically relevant range for that observable. 

The following has been added compared to the original version of the paper: The terms of 
the form dcr^^yj^Q in the five-parton and four-parton contributions to the N"^ and colour 
factors were only implemented in this revised version. They lead to changes in the numerical 
values of the NNLO coefficients which are most pronounced in the approach to the two- 
jet region. In a recent work, Becher and Schwartz [A] have computed the logarithmically 
enhanced terms which dominate the thrust distribution in the two-jet region using soft- 
collinear effective theory. They identified a disagreement with our original numerical results 
for the thrust distribution in the two-jet region for these two colour factors. Our new 
results are displayed in Figures 9, 10 and 11, and are now in full agreement with the results 
obtained in [A]. Our numbers also agree with the results obtained in the implementation 
of [B]. 

In the genuine three-jet region, which is relevant for precision phenomenology, the 
changes have a minor numerical impact. The corrections to the NNLO and colour 
factors also affect all other event shape distributions [C] in a similar manner; minor nu- 
merical effects in the three-jet region, but more significant effects in the two-jet region. 

15. Conclusions and Outlook 

In this paper, we provide a detailed description of the calculation of NNLO QCD corrections 
to three-jet production and related event shapes in electron-positron annihilation. At this 
order, three-parton, four-parton and five-parton subprocesses contribute. The three-parton 
and four-parton subprocesses contain explicit infrared singularities from loop corrections. 
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Four-parton and five-parton subprocesses contain singularities which only become explicit 
after integrating the contributions over the phase space relevant to the three-jet final states. 
Those singularities arise when one or two partons become unresolved (coUinear or soft). 
For an infrared-safe observable, adding together all infrared singularities, one observes a 
complete cancellation, resulting in an infrared-finite result. 

To implement the four-parton and five-parton processes in a numerical programme, 
one has to devise a procedure for extracting the implicit infrared singularities from them. 
We extract these singularities using subtraction terms which numerically subtract all in- 
frared singularities from the five-parton and four-parton channels. The subtraction terms 
are then integrated analytically and combined with the three-parton channel, where they 
cancel all explicit infrared poles. The subtraction terms are derived using the antenna sub- 
traction method, which is based on antenna functions encapsulating all unresolved partonic 
radiation emitted from a pair of hard radiator partons. 

Three-jet production at NNLO receives contributions from seven different colour fac- 
tors. Among these, only the six colour factors of non-singlet configurations require subtrac- 
tion, while the singlet colour factor is separately finite in all three partonic channels. We 
describe the construction of the antenna subtraction terms for the six non-singlet colour 
factors in detail, and demonstrate the cancellation of infrared poles. 

All partonic channels have been implemented in a parton-lcvcl event generator pro- 
gramme EERAD3, which can be used to compute any infrared-safe observable related to 
three-jet final states in e+e~ annihilation. We devised various tests of the implementation, 
demonstrating in particular the correct numerical cancellation between matrix elements 
and subtraction terms and the independence on phase space restrictions deep inside the 
subtraction regions. 

We observe that the largest part of the NNLO correction is contained in the two 
leading colour factors A^^ and Np N. The remaining four colour factors yield corrections 
at or below the ten per cent level. 

First phenomenological results on the thrust distribution at NNLO were obtained 
already in an earlier paper [33]. In the thrust distribution, the NNLO corrections amount 
to about 15% of the total result. They are lower in magnitude than the NLO corrections, 
indicating the perturbative stability of this observable. Inclusion of the NNLO corrections 
considerably reduces the dependence of the result on the renormalisation scale. 

At LEP, a wide variety of QCD event shapes was measured to high precision [70]. 
The accurate extraction of the strong coupling constant ag from these data sets was up 
to now limited by the theoretical uncertainty inherent to the available NLO calculations. 
We expect that our new NNLO results will improve this situation considerably. Numerical 
studies of other event shape variables and of the three-jet rate are ongoing, and will be 
reported elsewhere. 

On approaching the two-jet limit ((1 — T) — in the thrust distribution), one observes 
that the perturbative fixed-order expansion starts to break down due to the emergence of 
large logarithmic corrections at all orders in perturbation theory. By matching NLO calcu- 
lations and NLL resummation of these logarithmic corrections [11], a reliable description of 
distributions over the full kinematic range could be accomplished. Such a matching is also 
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possible at NNLO, where it requires the derivation of a number of new matching constants 
for each event shape. 

The subtraction terms derived here for e+e~ 3 jets at NNLO can be transcribed to 
the crossed reactions ep — (2 + 1) jets and — F + 1 jet at NNLO without much mod- 
ification. In these cases, which involve partons in the initial state [41], the same antenna 
functions are used with different antenna phase spaces. To accomplish the above-mentioned 
NNLO calculations therefore still requires the analytical integration of the relevant antenna 
functions over the phase spaces relevant to their initial-state kinematics, which appears fea- 
sible with present technology. NNLO calculations of other exclusive observables at hadron 
colliders, such as pp ^ 2 jets, could also be carried out using the antenna subtraction 
method by constructing their subtraction terms along the lines described here. 
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